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ABSTRACT 

We use 3 years of data from the Swift/BAT survey to select a complete sam- 
ple of X-ray blazars above 15keV. This sample comprises 26 Flat-Spectrum Ra- 
dio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 
0.03<z<4.0. We use this sample to determine, for the first time in the 15-55 keV 
band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs 
detected by BAT, the population of blazars shows strong positive evolution. This 
evolution is comparable to the evolution of luminous optical QSOs and luminous 
X-ray selected AGNs. We also find evidence for an epoch-dependence of the evo- 
lution as determined previously for radio-quiet AGNs. We interpret both these 
findings as a strong link between accretion and jet activity. In our sample, the 
FSRQs evolve strongly, while our best-fit shows that BL Lacs might not evolve at 
all. The blazar population accounts for 10-20 % (depending on the evolution of 
the BL Lacs) of the Cosmic X-ray background (CXB) in the 15-55 keV band. We 
find that FSRQs can explain the entire CXB emission for energies above 500 keV 
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solving the mystery of the generation of the MeV background. The evolution of 
luminous FSRQs shows a peak in redshift (zc=4.3±0.5) which is larger than the 
one observed in QSOs and X-ray selected AGNs. We argue that FSRQs can be 
used as tracers of massive elliptical galaxies in the early Universe. 

Subject headings: cosmology: observations - diffuse radiation - galaxies: active 
X-rays: diffuse background - surveys - galaxies: jets 



Introduction 



Blazars constitute the most extreme class of active galactic nuclei (AGNs). Their broad- 
band and highly variable eraission is due to a relativistic jet pointing close to our line of 



sight (e.g. iBlandford fc Reeslll978l ). In the framework of the AGN unified model, which as- 



cribe s the observed features of AGNs to orientation effects (jAntonucci 1993 : Urry fc Padovani 
19951 ). the properties of misaligned blazars are consistent with those of radio galaxies. In- 
deed, the two blazar sub-populations, BL Lacertae (BL Lac) objects and flat spectrum radio 
quasars (FSRQs), are thoug ht to be the beamed counterparts of low- and high- luminosity 
radio galaxies, respectively (IWall &: JacksonI 119971 : IWillott et al.ll200ll ). Both classes of ob- 
jects are normally fo und only in the nuclei of giant elliptical galaxies (e.g. lFalomo et al.ll2000l : 
O'Dowd et aLll2nn2h . 



Blazars have been extensively stud i ed at radio (jDunlop &: Peacock 



20051). soft X-ray (IGiommi fc Padovani 



1994 



Rector et al 



Caccianiga et al.ll2002[TBeckmann et al.lboosl : iPadovani et al.l 



2000 



1990 



Wall et al. 



Wolter fc Celotti 



20071 ) and GeV energies (IHartman et al 



2001 



19991 ). It seems consolidated that FSRQs evolve positively (i.e. there were more blazars in 



the past,lDunlop fc Peacock 



Padovani et al. 



200 



1 



Wall 



19901) up to a redshift cut-off which depends on luminosity (e.g. 



20081 ). In this respect FSRQs evolve similarly to t he population of 



X-ra y selected, radio-quiet, AGNs (lUeda et al.ll2003l : iHasinger et al.ll2005l : iLa Franca et al. 



20051 ). On the other hand, the evolution of BL Lac obje cts remains a matter of debate, since 



they were found to evolve negatively in a fe w cases (e.g. 



20031 ) and not evolving at all in other ones (jCaccianiga et al 



Rector et a. 



20021: 



200 



Beckmann et al. 



Padovani et al 



200 



3). 



Deriving the luminosity function of a class of objects allows to understand the properties 
of the parent population and to estimate the diffuse (unresolved) background produced by 
the entire class. Despite all the previous studies, the lack of a sensitive all-sky hard X- 
ray survey has prevented, so far, to gather a sizable sample of blazars and to study their 
evolution in the >10keV band. For these reasons, the contribution of blazars to the X-ray 
background (above > lOkeV) has never been quantified. The aim of this study is to address 
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all these questions using data from the Swift/BAT instrument. 



1.1. The high-energy Background 



Radio-quiet AGNs are more abundant than blaz ars and have beeii shown to be the maior 



constituent of th e Cosmic X-ray Background (CXB, lUeda et al.ll2003l : iTreister fc Urryll2005 



Gilli et al.ll2007l ). This consolidates the idea that the CXB emission is the result of accretion 
onto super-massive black holes. More precisely, the X-ray emission of AGNs is due to 
Compton up-scatter of UV photons (generated in the inner part of the accretion disk) by 
high-energy electrons which populate a region above the disk cor amonly ref e rred t o as corona. 
This process, known as Comptonization, was first proposed by IZdziarskil (119861 ). The bulk 
of the electron population present in the corona is expected to be thermal and this naturally 
produces a cut-off in the spectrum of AGNs. Th e detection of the AGN cut-offs, detected at 
energies between 50keV and 400 keV by OSSE jMadeiski et al.lll995l : Izdziarski et"aDl200ol 



and references therein) and th e non-detection of Seyfert-like AGNs by EGRET (ILin et al. 
19931 : iDermer fc Gehrelslll995l ) confirms this interpretation. 



Ueda et al. 


2003: 


Gilh et al. 


2007) 



the cut-off combined with the cosmic evolution of AGNs implies that the the contribution 
of radio-quiet AGNs to the CXB emission above ~200keV is negligible. Thus, the high- 
energy CXB emission in the 200-10000 keV energy range remains currently unexplained. A 
few candidates have been proposed to explain this background. O ne is the 7-ray ernission 



originated from nuclear decays froni Type la supernovae (SNe la; IClayton fc WardI Il975 



Zdziarskilll996l : IWatanabe et al.lll999l ). However, on the basis of measurements of the cosmic 
SN la rates, recent studies showed that the background flux e xpected from SNs la is about an 



order of magnitude lower than the observed CXB emission (lAhn et al.ll2005l : IStrigari et al. 



20051 ). Annihilation of dark matter particles has also been discuss ed, but no viable light 



(with "MeV" mass) dark matter particle candidate has been found (lAhn fc Komatsul 12005 
and references therein). 



Very recently llnoue et al.l (120081 ) discussed the possibility that the hot corona may con- 
tain a small fraction (relative to the whole population) of non-thermal electrons. These 
electrons might be powered by magneti c reconnections in a similar way as it happens dur- 
ing solar flares (e.g. IShibata et al.lll995l ). In this framework a faint non-thermal component 
present in millions of AGNs might explain the observed background. However, due to the 
lack of sensitive instruments surveying the MeV sky this hypothetical non-thermal emission 
of AGNs has never been detected. 
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Blazars, whose emission extend from the Radio to the TeV band, are certainly con- 
tributing, despite their relative low space density, to the high-energy background (both in 
the X-ray and 7-ray energy band). In particular, an important role is certainly played by 
the so called "MeV blazars" whose Inverse Compton (IC) pea k is located in the MeV band 
(IBloemen et al.lll995l : ISikora et al.ll2002l : ISambruna et al.ll2006l ). An attempt to quantify the 
contri bution of blazars to the high-energy background has been performed by iGiommi et al. 
(I2OO6I ) using a multi-frequency selected sample of blazars. Despite the uncertainties related 
to the extrapolation from the microwave to the hard X-ray energy band, they conclude that 
blazars can explain ~10 % of the CXB emission in the 2-10 keV energy band and possibly 
100% of the background above 500 keV. 

In this work, we use a complete sample of blazars detected in the Swift/BAT survey 
to derive, for the first time at these energies, the X-ray luminosity function (XLF) and the 
cosmic evolution of blazars and to assess their contribution to the high-energy background. 
The paper is organized as follows. In § |2] we describe how the blazar sample was selected 
among the BAT extragalactic sources and discuss its incompleteness. In § [3] we introduce 
the Maximum Likelihood method which is used to determine the blazar evolution. The 
luminosity function of the BAT blazars is derived in § H] and in § E] it is used to quantify the 
contribution of blazars to the diffuse background. We discuss the results of our analysis in 
§ [3 Throughout this paper we assume a standard concordance cosmology (Ho=70km s~^ 
Mpc-\ fiM=l-f^A=0.3). 



2. The Swift/BAT sample 



The Burst Aler t Telescope (BAT; iBarthelmy et al.l |2005| ) onboard the Swift satellite 
( IGehrels et al.l 120041 ). represents a major improvement in sensitivity for imaging of the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° 
partially coded) aperture sensitive in the 15-200 keV domain. The main goal of BAT is to 
locate Gamma- Ray Bursts (GRBs). While chasing new GRBs, BAT surveys the hard X-ray 
sky with an unprecedented sensitivity. Thanks to its wide FOV and its pointing strategy, 
BAT monitors continuously up to 80% of the sky every day. Thanks to this quasi-random 
pointing strategy and the large FOV, BAT exposure is uniform on the whole sky. 

^Results of the BAT survey jMarkwardt et all boosl . Aiello et all boOSai iTueller et al" 

2OO9I ) show that BAT reaches a sensitivity of ~1 mCraqj in 1 Ms of exposure. Given its 
sensitivity and the large exposure already accumulated in the whole sky, BAT poses itself as 



^ImCrab in the 15-55 keV band corresponds to 1.27x10 ^'^erg cm ^ s 
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an excellent instrument for studying populations whose emission is faint in hard X-rays. 

For the analysis presented here we used Swift/BAT survey observations performed be- 
tween March 2005 and Mar ch 2008. Dat a screen ing and processing was performed according 



to the recipes presented in lAjello et al.l (j2008al ). The chosen energy interval is 15-55 keV. 



The lower limit is dictated by the energy threshold of the detectors. The upper limit was 
chosen as to a void the presence o f strong background lines which could worsen the overall 



sensitivity (see lAjello et al.ll2008bl . for details about the BAT background). The all-sky im- 
age is obtained as the weighted average of all the shorter observations. The average exposure 
time in our image is 4.3 Ms, being 2.0 Ms and 6.8 Ms the minimum and maximum exposure 
times respectively. The final image shows a Gaussian normal noise and we identified source 
candidates as excesses above the 5 a level. All the candidates are then fit with the BAT 
point spread function (using the standard BAT tool batcelldetect) to derive the best source 
position. Moreover, in order to avoid problems related to source confusion and sample in- 
completeness, we considered only sources at high (|b|>15°) Galactic latitude. This analysis 
is based on mean source fluxes determined over the 3 year period spanned by the survey. 
Our high-latitude sample comprises 305 sources. Of these 40 are Galactic sou rces (mainly 



X-ray binaries) and 6 are galaxy clusters (already comprised in the sample of lAjello et al. 



( 120091 )). The remainder are 247 extragalactic sources and 12 unidentified objects. The exact 



composition of the sample is reported in Tab. [H 

Being the BAT survey not a flux-limited survey, but rather a significance-limited one, 
it is important to address how the survey flux limit changes over the sky area. This is often 
referred to as sky coverage, that is the distribution of the survey's area as a function of 
limiting flux. Its knowledge is very important when perfor ming popula t ion stu dies as the 



ones described in the next sections. The reader is referred to lAjello et al.l (l2008al ) for how to 
derive the sky coverage which as a function of the minimum detectable flux Fmin is defined 
as the sum of the area covered to fluxes fi < Fmin'- 



N 



^ ( ^ Fmin ) ^ ^ -^i 1 fi'^ Pmin ( 1 ) 

i 

where N is the number of image pixels and Ai is the area associated to each of them. 
A visual representation of the sky coverage is reported in Fig. [T] which shows clearly the 
good sensitivity of BAT which reaches, in our analysis (15-55 keV), a limiting sensitivity of 
~ 0.6mCrab (7.3x 10~^^ erg cm"^ s~^). 
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Table 1. Composition of the BAT high-latitude sample (|b| >15° and S/N>5). 



CLASS 


# ol")jects 


Total 


305 


Seyferts 


199 


Blazars 


38 


Galaxies 


4a 


Galaxy Clusters 


6 


Radio Galaxies 


6 


Galactic Sources 


40*^ 


Unidentified 


12 



'^These objects are candidate 
radio-quiet AGNs (Seyferts), 
identified by means of a 2- 
10 keV follow-up observation, 
for which an optical spectrum, 
and thus redshift, is not yet 
available. 

^It includes all objects of 
Galactic nature (i.e. pulsars. 
X-ray binaries, etc.) 
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Fig. 1. — Sky coverage of the BAT survey as a function of minimum detectable flux for 
S/N> 5 a and |b| >15° in the 15-55 keV band. 

2.1. The Blazar sample 



Selection of blazars in the BAT energy band is not problematic. Thanks to the deep, 
although sparse, coverage of the sky in soft (0.1-2.5 keV) and medium (2-10 keV) X-rays, 
most of the BAT sources are already we ll stud ied objects. The location accuracy of BAT was 
recently characterized by iTueller et al.l (120091 ) using the all-sky sample of sources detected 
in 22 months of observations. They find that the 96 % error radius for a 5 a, 10 a, and 
20(7 source is respectively 7.5' , 3.34' and 1.59' . For the identification of most of the BAT 
sources, we used the R QSAT All-Sky Su rvey Bright Source C atalogue (IVoges et al.lll999l) . 
the Third IBI S Cata log feird et al.l l2007h. the BAT catalogs of iMarkwardt et al.l toO^ and 
Tueller et aP J2OO8I . hoO^ as well as SIMBAD and NED. We remark that all the identi- 
fications reported in the aforementioned catalogs are based on a thoro ugh study, both in 



the optical and X- rays, of the properties of the BAT counterparts (e.g. lAjello et al.ll2008c : 
Winter et al.ll2008l ). Some 50 BAT sources, for which the above catalogs did not provide 
an identification, have a publicly available Swift/XBT observation. These observations pro- 
vided a secure identification of the BAT objects and the results of all these follow-ups will 
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be discussed elsewh ere (see Tab. [T] for a br eakdown of the total sample) . As shown in 
Aiello et all J2008c[ ) and Iwinter et al.l boO^ . most of the BAT AGNs show complex X-ray 
spectra (with the presence of iron line, Compton reflection, and soft excesses) and often very 
large absorption (Nh ~10^^ atoms cm~^) which securely associates the BAT object with a 
non-blazar AGN (e.g. Seyfert). 

In order to securely identi fy and classify BAT sources as blazars, we relied mainly on 
the blazar catalog (BZCAT) of iMassaro et al.l (120091 ) which contains only bona-fide blazars. 
These are sources which show a dominant, broad-band, non-thermal component associated 
with a (relativistic) jet, a variable nuclear activity and are detected at least in radio, optical 
and X-rays. Given the fact that the BZCAT do es not cover the en tire sky, we also used 
the CRATES catalog of FSRQs (see also below; iHealey et al.l 120071 ) and when necessary 
individual source publications. The sample of BAT blazars is reported in Tab. [2] along with 
the main properties of the sources (e.g. fluxes, signal-to-noise ratios, etc.), while the sample of 
Seyferts and radio galaxies detected by BAT is reported for reference in the Appendix |Al For 
this analysis, BL Lacs are identified as objects in which the equivalent width of the strongest 
emission line is less than 5 A and the optic al spectrum shows a Ca II H/K break <0.4 (e.g. 
Urry fc Padovanilll995l : iMarcha et al.lll996l ). The blazar sample comprises 26 FSRQs and 12 
BL Lac objects and the full details of the classification (including the references for it) are 
given in Tab. [2l Among the BL Lac objects, 9 are of the high-frequency-peaked (HBL) type 
while the rest are of the low-frequency peaked (LBL) type. Figure [2] (both panels) shows that 
the BL Lacs detected by BAT have a softer spectrum with respect to the FSRQ population. 
The average photon index of the BL Lac objects, in the BAT band is 2.5±0.5 while the one 
for FSRQs is 1.6±0.3. This is however not surprising, indeed it confirms the expectation that 
BAT samples the synchrotron component in BL L acs and the IC component in FSRQs (e.g. 
see the BeppoSAX results of iDonato et al.ll2005l ). Figure [3] shows the luminosity- redshift 
plane for the BAT blazars in comparison with the Seyfert-like AGNs detected by BAT. 
There are a few things that can be noted. First, the whole BAT sample of AGNs (Seyferts 
and blazars) spans almost 4 decades in redshift and 8 in luminosity. There are 16 objects 
detected at redshift larger than one and they are all FSRQs. On the other hand, no blazars 
are detected at low luminosities and low redshift. This will be discussed in the next sections. 



We cross-correlated the BAT blazar sample with the Third EGRET catalog (IHartman et al 



19991 ) using a two degree search radius around the BAT positions to cope with the large 
error radius of the EGRET positions. We further made sure, by searching in the litera- 
ture, that the EGRET source , if any, is securely associated with the BAT counterpart (e.g. 
Sowards-Emmerd et al.l |2004| ) . We found that only 12 BAT sources are included in the 
EGRET catalog being these 3 BL Lac and 9 FSRQ objects (see Tab. [2]). Moreover, we also 
checked how many BAT blazars were also detected by Fermi-LAT in the first three months 
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Redshift 



2 2.5 
Photon Index 



Fig. 2. — Left Panel: Photon index versus redshift for the blazars in the BAT sample. 
The low redshift sources, mainly BL Lacs (red squares), have a soft spectrum (photon index 
of ~2.5, while the high redshift ones (mainly FSRQs, blue circles) have a hard spectrum. 
Right Panel: Photon index distribution for FSRQs (solid line) and BL Lacs (dashed line). 



of operations (lAbddl2009bl ). Again we found that only 12 sources (9 FSRQs and 3 BL Lacs) 
are in common between the two samples (see Tab. [2]). This implies that EGKET/ Fermi 
and BAT do not sample exactly the same blazar population. In particular, the blazars BAT 
detects at the highest redshifts' are not (yet) present in the EGKET /Fermi catalogs. This 
would mean that the high-redshift BAT blazars are faint in the GeV band. This hypothesis 
will be discussed in the next sections. 

We believe the incompleteness of our blazar catalog is negligible. The fraction of uniden- 
tified sources in the total extragalactic sample is <5 %. The most simple assumption is that 
these sources are distributed in source classes exactly like the identified sources. However, 
this does not take into consideration other properties like rad io emission. Indeed, no blazars 
are known to be radio-quiet sources (jWolter fc Celottill200lh . On this basis we found that 
none of the unidentified sources has a radio-loud object within 8' or is present in the blazar 
catalogs used above. In addition we also looked for spatial coi ncidences, within 8 ^ , be- 
tween the BAT unidentified objects and the CRATES catalog (IHealey et al.l 120071 ) . The 
CRATES catalog contains all known FSRQs object above a 4.8-GHz fiux of 65mJy and 
for |b|>10 degrees. Again we found that none of the BAT unidentified objects is associated 
with a CRATES radio source. We, thus, confidently believe that the incompleteness of the 
BAT blazar sample is negligible (i.e. no new blazars are hiding among the unidentified BAT 
sources) . 



Table 2. Blazar sample 



SWIFT NAME 


R.A. 
(J2000) 


Dccl. 

(J2000) 


Flux 

(10~^^ erg cm"^-^ s~^) 


S/N 


ID 


Type 


Rcdshift 


Photon index 


3EG^ ? 


LAT ? 


Rcf. 




JOOIO. 4+1056 


2.617 


10.935 


1.88+0.21 


9.2 


QSO B0007+107 


FSRQ 


0.09 


2.08 + 0.40 






B 




J0018. 8+8137 


4.713 


81.624 


1.18+0.19 


6.2 


S5 0014+81 


FSRQ 


3.36 


1.93+0.55 






B 




J0123. 0+3421 


20.752 


34.351 


1.27+0.21 


5.9 


lES 0120+340 


BLLAC (HBL) 


0.27 


3.05+0.60 






B 




J0207. 3+2929 


31.848 


29.500 


1.25±0.22 


5.6 


4C 29.06 


BLLAC (LBL) 


0.11 


2.52+0.65 






V 




J0225. 2+1849 


36.302 


18.823 


1 .42 + 0.22 


6.5 


RBS 0315 


FSRQ 


2.69 


1.48+0.45 






R 




J0233. 1+2017 


38.292 


20.295 


1.65+0.22 


7.4 


lES 0229+200 


BLLAC (HBL) 


0. 14 


2.33+0.50 






B 




J0313. 0-7645 


48.250 


-76.750 


0.97+0.12 


5.1 


PKS 0312-77 


FSRQ 


0.22 


1.66+0.97 






C 




J0336. 5+3219 


54.136 


32.325 


1.83+0.24 


7.6 


4C 32.14 


FSRQ 


1.26 


2.03±0.50 






B 




J0349. 7-1157 


57.449 


-11.951 


1 .32+0.21 


6.3 


lES 0347-121 


BLLAC (HBL) 


0.19 


2.36±0.35 






S 




J0353. 1-6829 


58.287 


-68.490 


1.34+0.18 


7.5 


IGR J03532-6829 


BLLAC (HBL) 


0.09 


3.33d=0.81 






I 




J0523. 0-3626 


80.755 


-36.447 


1.61±0.18 


9.0 


PKS 0521-365 


BLLAC (LBL) 


0.06 


1.92+0.38 


y 




F 




J0525. 5-4557 


81.398 


-45.951 


1.06+0.17 


6.4 


PKS 0524-460 


FSRQ 


1.48 


1.31+0.38 






K 




J0539. 9-2838 


84.999 


-28.650 


1.27+0.20 


6.2 


PKS 0537-286 


FSRQ 


3.10 


1.56+0.30 






C 




J0550. 8-3217 


87.716 


-32.289 


2.08+0.20 


10.5 


PKS 0548-322 


BLLAC (HBL) 


0.07 


1.88+0.28 






S 




J0635.9-7515 


99.000 


-75.250 


0.94+0.18 


5.1 


PKS 0637-752 


FSRQ 


0.64 


1.87+0.57 






c 




J0746. 5+2550 


116.648 


25.848 


1.49+0.25 


5.9 


SDSS J074625. 87+254902. 2 


FSRQ 


2.98 


1.08±0.38 






B 




J0805. 3+6148 


121.349 


61.800 


0.96+0.19 


5.1 


GB6 J0805+6144 


FSRQ 


3.03 


1.58±0.62 






B 




J0841.4+7054 


130.363 


70.915 


2.85+0.18 


15.9 


4C +71.07 


FSRQ 


2.17 


1.47+0. 14 


y 




B 




J1104. 5+3812 


166. 126 


38.210 


12.16+0.17 


80.5 


Mrk 421 


BLLAC (HBL) 


0.03 


2.64+0.10 


y 


y 


B 




J1130. 0-1448 


172.512 


-14.815 


2.17+0.26 


8.3 


PKS 1127-145 


FSRQ 


1.19 


1.99=1=0.39 


y 


y 


T 




J1213. 2+3238 


183.300 


32.648 


0.90+0.17 


5.4 


B2 1210+33 


FSRQ 


2.50 


1.40=t0.30 






B 


h 

c 


J1224. 9+2118 


186.249 


21.300 


0.95+0.18 


5.3 


QSO B1222+216 


FSRQ 


0.43 


1.53±0.54 


y 




B 


J1229. 1+0202 


187.283 


2.047 


18.31+0.21 


87.0 


3C 273 


FSRQ 


0.16 


1.70+0.07 


y 


y 


B 




J1256. 1-0547 


194.048 


-5.800 


1.40+0.23 


6.1 


3C 279 


FSRQ 


0.54 


1.61+0.35 


y 


y 


C 




J1428. 8+4240 


217.208 


42.667 


1.40+0.16 


8.5 


H 1426+428 


BLLAC (HBL) 


0.13 


2.70+0.50 






B 




J1513. 1-0903 


228.275 


-9.061 


2.49+0.34 


7.3 


PKS 1510-089 


FSRQ 


0.36 


1.57+0.29 


y 


y 


B 




J1654. 1+3945 


253.542 


39.763 


3.46+0.21 


16.7 


Mrk 501 


BLLAC (HBL) 


0.03 


2.47+0.22 




y 


B 




J1959. 7+6509 


299.940 


65.159 


2.29+0.20 


11.5 


lES 1959+650 


BLLAC (HBL) 


0.05 


2.51+0.29 


y 


y 
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Fig. 3. — Luminosity- redshift plane for the BAT blazars (empty circles) and BAT Seyfert-like 
AGNs (filled gray circles). 



We first test t 



3. The evolution of Blazars 

le ev olution of the BAT blazars by applying the V^/Vmax method pro- 



posed by ISchmidtl (Il968l ). For a non-evolving source population, V/Vmax is expected to be 
uniformly distributed between and 1. Thus, in case of no evolution the average V/Vuax 
is 0.5. The error on the average value can be computed as cr = 1/(12A^)^/^. A value of 
the <V/Vmax> significantly different from 0.5 indicates positive (if > 0.5) or negative evo- 
lution (otherwise). Computing the <V^/Vmax> for the whole sample of blazars we obtain 
0.666±0.045. This indicates that the blazars population detected by BAT evolves positively 
at >3 cr. This means that the luminosity or density of the blazars is increasing with redshift. 
Moreover, we computed the <V^/Vmax> for FSRQs and BL Lacs separately in order to test 
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if the two sub-populations evolve differently. We find a value of <V /Vmax> of 0.728±0.056 
and 0.576±0.083 for FSRQs and BL Lacs, respectively. Thus, this preliminary analysis shows 
that while FSRQs are evolving strongly BL Lacs show only mild evolution (compatible at 
~ 2cr with no evolution). As a control sample we us ed the sample of the Seyfert-like AGNs 



detected by BAT. Since this population is truly local (jAjello et al.ll2008d : iTueller et al.ll2008l ) 
we expect the <V /VMAyi> to return a test value of 0.5. As expected, we obtain 0.509±0.021. 
These results are summarized in Table [31 

Another test of cosmological evolution, although less powerful than the l^/HiAx method, 
is the log iV-log S test which is based on the distribution of source counts above a given flux. 
This distribution is generally expressed as: A^(> S)=A. where S is the source flux. In this 
test, a non-evolving population shows a distribution which is consistent with an Euclidean 
distribution with /5 = 3/2. If the source population is positively evolving then a /3 > 3/2 is 
expected. The results of the log A^-log S test are in excellent agreement with those of the 
<V^/Vmax> method (see Table [3]). They confirm that while the Seyfert population is not 
evolving, the blazar populations (and in particular the FSRQs) are evolving positively. Also 
the normalization of the log A^-log S distributions (e.g. A parameter in Table [3l) highlights 
that, at the current flux limit of BAT, blazars are 5 times rarer than normal Seyfert galaxies. 
Thus, it is only thanks to the uniformly deep all-sky exposure that BAT gathered a sizable 
sample of them. 



3.1. The evolving Luminosity Function 

The differential luminosity function of a population of objects is defined as the number 
of objects per unit comoving volume and per unit luminosity interval: 

<fN , , dV , , dV 
where dV/ dz is the comoving volume element per unit redshift and unit solid angle (see 



e.g. lHoggill999l ). The present day XLF can be obtained for z=0 as $(I/x, z = 0). We model 
the present day XLF with two commonly used functions. A simple power of the form: 
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and a double power-law of the form (see e.g. lUeda et al.ll2003l : iHasinger et al.ll2005l ): 



dN 



A 



dLx /n(10)Lx 



71 



+ 



72 



(4) 



Where in Eq. [3] we set = 10^^ erg s ^ while in Eq. |4] is allowed to vary. 



3.1.1. The luminosity and density evolutions 

The simplest scenarios of evolution are pure luminosity (PLE) and pure density evolution 
(PDE). In the PLE case the XLF becomes: 



^Lx{z),z) = ^{Lx/e{z),z = Q) 



(5) 



while in the PDE case: 



^{Lx,z) = ^{Lx,z = Q) X e{z) 



(6) 



and the evolution is parametrized using the common power-law evolutionary factor: 

e{z) = (1 + zf^^^ (7) 



It is thus clear that in the PLE case the typical blazar luminosity is changing with 
redshift while in the PDE case only their densities are changi ng with red s hift. We note that 
this formulation of the evolutionary factor (suggested first by IWall et al.l (120081 )) reduces for 
7 = to the standard power law (1 -|- zY . For clarity, we call PDE (PLE) or modified- 
PDE/PLE (or MPDE, MPLE) the case in which 7 has been fixed to zero or was allowed to 
vary. 



3.2. Maximum Likelihood Analysis 



A classical approach to derive the XLF is based on the 1/Vmax method of ISchmidt 



(Il968l ). However, this method is known to introduce a bias if there is significant evolution 
within bins of redshift. Moreover, considering the small number of objects in our sample, 
binning would result in a loss of information. We t hus decided to apply the Maximum 
Likelihood method using the formalism introduced by [Marshall et al.l (119831 ). 
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In this method, the luminosity-redshift plane is parsed into extremely small intervals of 
size dLxdz. In each element we compute the expected number of blazars with luminosity 
Lx and redshift z: 



dV 

\{Lx, z)dLxdz = p{Lx, V)n{Lx, z)— dLxdz (8) 

dz 

where VL{Lxiz) is the sky coverage of the survey (see §[2] for details). The sampling of 
the luminosity-redshift plane is sufficiently fine that in each dLxdz element the number of 
observed blazars is either 1 or 0. In this sparse sampling limit we can define a likelihood 
function based on joint Poisson probabilities where the Poisson model is: 



f{x : m) 



x\ 



(9) 



where /x is the expected number of blazars. If x = 1, then the function is /ie ^ and if x = 0, 
it is e~^. In this case the likelihood function can be written as: 



L = Y[\{Lx,i,Zi)dLxdze-^^^'''"'''^'^^'"^' x JJ g^^^^^'^'"^)'^^^'^^ 



(10) 



This is the combined probability of observing one blazar at each element {Lx,i, Zi) popu- 



lated by one BAT blazar and observing zero blazars everywhere else {Lx. 



31 ^3 1 



Transforming 



to the standard expression S = —2ln L and dropping terms which are not model dependent, 
we obtain: 



S = -2J2w^ Hp{Lx,i,Zi)]+2 



X{Lx,i,Zi)dLxdz 



(11) 



where following iBorgani et al.l (120011 ) we have introduce a weighting term Wi which takes 
into account the uncertainties in the luminosity of each single blazar. In this way, each blazar 
instead of being a point in the Lx,^-plane is smoothed in the Lx-direction according to a 
Gaussian distribution with a width set by the 1 a luminosity error. Thus, a weight is assigned 
to each element in the luminosity-redshift plane based on the fractional contributions of all 
blazars in the same redshift interval: 



k 



1 



27re 



exp 



{Lx,k — Lx^iY 



2e2 



LX,k 



dL 



X 



(12) 



where the summation k is over the blazars with a redshift between Zi — dz/2 and Zi + 
dz/2. The limits of integration of Eq. [11], unless otherwise stated, are: Lx,min=^0'^^ erg s~^. 
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-Z^x,max=10^° erg s~^, Zmin=0 sjid Zmax=G- While, the results are independent of the upper 
limits of the integration, as far as they are chosen to be large, this is not the case for the 
value of Lx,min if the local XLF is modeled as a single power law. In this case, the lower 
limit of integration needs to be set to the minimum observed blazar luminosity (6x10^^ erg 
s~^ and 2x10^^ erg s~^ for BL Lacs and FSRQs respectively). 



The best fit parameters are determined by minimizingj S and their associated 1 a errors 
are computed by varying the parameter of interest, while the others are allowed to float, until 
an increment of A5'=l is achieved. T his gives an estimate of the 68% confidence region for 



the parameter of interest (jAvnilll976l ). 



3.3. Alternative Maximum Likelihood Formulation 



Another way of posing the M aximum Likelihood problem is (following e.g. lChiang fc Mukherjee 



19981 : iNarumoto fc Totanill2006[ ): 



L = exp(-Ar,,p) Y[ ^Lx,^, Zi) (13) 



1=1 



where Nexp is the expected number of blazar detections: 

Ne.p = jdzj dLx^Lx,z). (14) 

In this case the function S (= —2lnL) is defined as: 

Nobs 

S=-2J2 HHLx^i, Zi)) - 2N ln{N,,p). (15) 

i 

We tested that we get exactly the same results if we use one or the other formulation 
of the Maximum Likelihood problem. Thus, the results that we present in the following 
sections are independent on the ML function chosen. 



3.4. Consistency Checks 

The Maximum Likelihood approach does not provide a goodness of fit test and this 
implies that other methods have to be used to understand if the fitted function is a good 



^The MINUIT minimization package, embedded in ROOT (root.cern.ch), has been used for this purpose. 
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representation of the data. A common procedure is to use the Kolmogorov-Smirnov (KS) 
test which is based on the maximum distance [D^s) between the (cumulative) distributions 
under comparison. This test computes the probabihty of observing a KS test statistics as 
large or larger than the observed one and it can be used to reject a model when too low. 
We apply a KS test to both the cumulative redshift and luminosity distributions of the BAT 
blazars and we reject XLF models which produce a KS probability <20 %. 

As a further test we check that the best-fit XLF reproduces well the observed source 
count distribution (also known as log - log S). The all-sky number of blazars with a flux 
stronger than S can be computed as: 

N{> S)=A7r dz— / dLxp{Lx, V{z)) (16) 

where Lx{z,s) is the luminosity of a blazar at redshift z whose flux is S. 



3.5. The Cosmic X— ray Background constraint 



It is almost certain that the bulk of the CXB emission (below 200 keV) , even if presently 



unresolved above 10 keV , is d ue to Seyfert-like AGNs (e.g. lUeda et al.ll2003l : iLa Franca et al. 



200,4 iTreister fc Urrvl l2005l : lOilh et all 120071 : [Silverman et al. 



200 



Even though one of 

our goal is to estimate, in the most robust way, the contribution of blazars to the CXB 
spectrum, as a first step the CXB emission can be used to reject invalid XLF models. 
Indeed the blaz ar contribution to t he CXB emission at X-ray energies is expected to be of 



the ~10 % (e.g. iGiommi et al.ll2006l ). A much larger fract ion would conflict with the present 



estirnates produced by population synthesis models (e.g. lUeda et al.ll2003l : ITreister &: Urry 



20051 : ICilli et al.l 120071 ) and can be used to rule out a given evolutionary model. Thus, for 



each best fit XLF model we compute the integrated background flux arising from the blazar 
population. This flux can be derived as: 

dV f^X.'rnax 



'CXB 



dz- 



dz 



dLxFx{Lx,z)p{L,,V{z)) 



(17) 



where the limits of integration are the same as in Eq. [TT]and Fx{Lx, z) is the flux of a 
source with luminosity Lx at redshift z. 

As a final note, we remark that we are not interpreting the CXB as a 'hard constraint' 
in the sense that the integrated blazar emission of Eq. [17] is not constrained to be ~10 %, 
but a model XLF which overproduces the entire CXB (e.g. producing more than 100 % of 
the CXB), in the 15-55 keV band, can be certainly ruled out. 
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Results 



4.1. The control sample: Local Seyfert galaxies 



Recently, iTueller et al.l (120081 ) computed the (non-evolving) XLF of local Seyfert galaxies 
using a sample of 88 objects detected by BAT in the first 9 months of operations. Given 
the small redshift range spanned (z<0.1), they did not test for the evolution of Seyferts in 
the local Universe. However, it is well established that the population of radio-quie t AGNs 
evolves in density and luminosity (e.g. iHasinger et al.ll2005l : [Silverman et al.ll2008l ). Thus, 



the sample of 199 Seyferts detected in this analysis represents a good test for the ML method 
introduced in § I3.2[ Here, we aim at deriving a parametric representation of the Seyfert XLF 
testing at the same time for their evolution in the local Universe. 



We model the local XLF (e.g. ^{Lx, z = 0)) with a double power-law model as in Eq. H] 
and fix the evolutionary term k, of the PLE, to zero. The best fit parameters, reported in 
Tab. S (model 1), are 71 = 0.80 ±0.08, 72 = 2.67 ±0.20 a nd L* = 6.1 ± 1 4 x 1 0^^ erg s'^ and 
are in very good agreement with the values reported by Tueller et al. (2008]f|. The error 
bars are generally smaller because the sample we use is larger and because the fit is done to 
the unbinned dataset. As shown by the redshift, luminosity and log A^-log 5* distributions 
(reported in Fig. [Hand [5]), this non-evolving XLF model is an highly acceptable description 
of the dataset (KS tests ~1, see Table H]). 

Allowing the XLF to evolve in luminosity (see model 2 in Tab. Hj) produces an equally 
good fit with an evolution parameter {k = 2.62 ± 1.18) which denotes positive evolution 
although constrained only at the ~2cr level. It is interesting to note that the evolution 
parameter is in good ag r eemen t wit h the values of 2.29± 0.09 and 2.7±0.2 found for the 
PLE case by lUeda et al.l (120031 ) and IHasinger et al.l (120051 ) respectively. However, since in 
our case the two models produce an equally good fit (see KS test values), the non-evolving 
XLF has to be preferred because of the lower number of free parameters. Thus, we believe 
that the evidences of the evolution of radio-quiet AGN in the local Universe are, with the 
current dataset, marginal. As a final proof, we built a non-para metric represen tation of the 
luminosity function of the Seyferts using the 1/Vmax method (ISchmidtlll968l ). In order to 



test for evolution we binned the dataset in two redshift bins containing approximately the 
same number of sources. This luminosity function is reported in Fig. [61 An evolving XLF 
would show a shift (in luminosity or density) from one redshift bin to the other. This is not 



■^The different value of Hq and energy band that Tueller et al. ( 2008h adopt produce the n et effect that 
the luminosities quoted here are directly comparable to those quoted by iTueller et al.l (j2008f ) without the 
need of a conversion factor (i.e. the conversion factor is 
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Fig. 4. — Redshift (left) and luminosity (right) distribution of the BAT Seyferts. Lo ng error 
bars consistent with zero are la upper limits in the case of observing zero events (IGehrels 
19851 ). For both cases, the dashed line represents a non-evolving XLF (model 1 in Tab Hj) 
convolved with the BAT sky coverage. 
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Fig. 5. — Cumulative log - log S of the BAT Seyferts. The dashed line is the prediction 
of the best-fit XLF (model 1 in Tab [4]). 



the case for the BAT Seyferts whose luminosity functions, derived in two different redshift 
bins, are the continuation one of each other. In the same Figure, the best-fit non-evolving 
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XLF (model 1) is also displayed and it is clear that this represents the data well. The fact 
that the evolution, if detected, is only marginal is not surprising, but consistent with the 
results of the V/Vm ax test repor t ed in T able p and the Euclidea n behavior of the log - 
log S as reported bv kiello et aP J2008ch and lTueller et al.l J2008I ). 
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Fig. 6. — Luminosity function of the Seyfert-like AGNs, derived using the 1/Vmax method, 
in two redshift bins (data points). The solid line is a non-evolving XLF double power-law 
model (model 1 in Tab. H]) as derived from the ML algorigthm. 



4.2. A single Blazar population 



As a first case, we start assuming that the local (present day) XLF ca n be adequately 
appro ximated by a simple power law as in Eq. [31 It was already shown by [Marshall et al. 
(119831 ) that in this case it is impossible to discriminate between density and luminosity 
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evolution. So we refer to density evolution, but we note that the two type of evolutions are 
formally equivalent in this case. 

In the simple power-law case (and PDE) we obtain a slope of the XLF of 72 = 2.67±0.13 
and an evolution parameter of 4.00±0.77. The error on the evolution parameter k confirms 
that the evolution is significantly detected in the BAT sample (see model 3 in Tab. H]). 
We note that the XLF slope is in very good agreement with the bright end slope of the 
Seyfert-like AGNs detected by BAT (see Tab. SJ. Figure [7] shows the redshift and luminosity 
distribution of the BAT blazars with superimposed the best fit XLF. The KS test shows that 
this is already an acceptable description of the data. 

We note, from left panel of Fig. [7l that the best-fit PDE XLF fails to describe the drop 
in blazar counts above z~4. This might be a sign of a possible cut-off in th e evolution. Thus, 



we decided to model the evolution factor as e{z) = (1 + 2;)^+'^^ as done by lWall et al.l (120081 ) 
(see Eq. [7]) calling this model a modified PDE (or MPDE). The best fit to the data shows 
that the value of 7 is constrained to be negative at the ~2.5(T level (7=-0.69±0.27). This 
shows that a cut-off in the evolution is needed in our data. Although, this model reproduces 
our data reasonably well, the integration of the XLF shows that it overpredicts the "total" 
CXB emission by a factor 3 (see Tab. H]). 




Fig. 7. — Redshift (left) and luminosity (right) distribution of the BAT blazars with super- 
imposed (dashed line) the prediction of the best-fit PDE XLF (single power-law case, see 
model 3 in Tab. Hj) convolved w ith the BAT s ky coverage. Error bars were computed taking 
into account the Poisson error (jGehrelslll985l ). Long error bars consistent with zero are la 
upper limits for the case of observing zero events in a given bin. 



Thus the interpretatio n of density evolution might not be the correct one. A similar 
effect was already noted by [Marshall et al.l (119831 ) who concluded for their sample of optical 
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QSOs that luminosity evolution was likely occurring. 

As already said, for a single power-law XLF there is no formal difference between density 
or luminosity evolution. The only difference is that in the integral of Eq. [TTl the integration 
limit Lx,min is evolving and can be expressed as Lx,min = Lxmin ^ ^i^)^ where L° 



X,min 



IS 



the present day luminosity cut-off (e.g. 6x 10^'^ erg s~^) and e{z) = (1 + z)^'^'^^. The best 
fit confirms that indeed the luminosity evolution is a better interpretation of the underlying 
evolution. Indeed, integrating the XLF we get that the blazar population, described by this 
MPLE function, accounts for 20 % of the CXB emission in the 15-55 keV band. We also note 
that the best-fit valu e of the evo l ution parameter k = 2.96 ± 0.46 is in ver y good agreement 



with what found by lUeda et al.l (120031 ) for X-ray selected AGNs and by IWolter fc Celotti 



tOQlh for X-ray selected FSRQs. 



All the XLF models described so far (models 3, 4 and 5 in Tab.Hj) become unacceptable 
if the limit on the minimum observed luminosity Lmin (see Eq. [TT!) is removed. Indeed, in 
this case (see results of the KS tests for model 6) the luminosity and redshift distributions are 
not reproduced correctly because the best-fit model predicts many blazars at low luminosity 
and low redshift which are not detected by BAT. Thus, a rather drastic change in the power- 
law behavior of the local XLF is required in order to reproduce the lack of low-luminosity 
objects. 

To test this scenario, we model the local XLF as a double power law (see Eq. H]) coupled 
to a MPLE model. In this model, we remove the constraint of a low luminosity cut-off and 
the fit is performed to the whole luminosity-redshift plane. This XLF model reproduces our 
data accurately (see model 7 in Tab. [Hand the distributions reported in Fig. [8]). Given the 
lack of low-luminosity objects, the faint-end slope 71 is, from our fit, required to be flat, but 
poorly constrained (-0.87±1.31). 

The slope of the bright-end part of the XLF 72 = 2.73 ± 0.37 is in good agreement with 
the same slope found for the Seyferts (see model 1). The likelihood ratio test can be used to 
assess whether a model produces a significant improvement over another one. The likelihood 
ratio test is the difference between the value of S (see Eq. [TT]) produced by different m odels. 



This value (AS*) is expected to be asymptotically distributed as the xi ( !Wilkslll938l ) where 
n is the difference between the degrees of freedom of the two models. The AS for model 7 
(double power law plus MPLE) with respect to model 4 (single power-law plus MPDE) is 
~10.1 which translates in a probability of 0.00150 that the improvement was obtained by 
chance. We also note that density evolution (with a double power law as a local XLF) 



*The chance probability was computed using the distribution for one degree of freedom. Indeed, the 
difference between model 7 and model 4 is given by the faint-end slope 71 which is allowed to vary while the 



-22- 



N 10 



10' 



lo- 



10"" 



10-^ 



10' 



1 10 



10 



10' 



10-* 



lo-* 



1. 



10-^ 10' 1 10 10^ 10^ 10' 10' 

[10* ergs'] 



Fig. 8. — Redshift (left) and luminosity (right) distribu tion of the BA T blazars. Error bars 
were computed taking into account the Poisson error (jGehrelslll985l ). For both cases, the 
dashed line represents the MPLE XLF (double power-law model, see model 7 in Tab. Hj) 
convolved with the BAT sky coverage. 



reproduces the data equivalently well (see model 8), but it is ruled out since it overpredicts 
the CXB emission by a factor > 5. Thus, we consider model 7 as the best representation of 
our data. Figure [3 shows the confidence contours for the best-fit parameters. 




Fig. 9. — Confidence contours (1,2, and 3 a) for the bright-end XLF slope (72), the evolution 
parameters {k and 7) and the break luminosity for the best-fit XLF model (model 7 in 
Tab.lH). 



The extrem e flattening of t he XL F at low luminosities can be the effect of beaming. 
As discussed by lUrry &: Shaferl (119841 ) relativistic beaming alters the observed luminosity 



cut-off luminosity (imposed in model 4) and the break luminosity L» (in model 7) represent essentially 

the same parameter. 
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function of bl azars producing a fla ttening at low luminosities. For common jet emission 



scenarios (see lUrry &: Shaferl Il984l . for details), the faint-end slope of the XLF should be 
~ 1.0. Given the absolute lack of BAT blazars populating the low-lumi nosity part of th e XLF , 
it is not surprising that the best-fit value of 71 is ~1.5 a away from the lUrry &: Shaferl (119841 ) 
prediction. On the other hand, relativistic beaming should not affect the bright end slope 
which should reflect the slope of the intrinsic luminosity function. It th us becomes inte r esting 
to compare the value of 72 derived here with other surveys. Recently ICara fc Lister! (120081 ) 
derived the intrinsic radio luminosity function of the Fanaroff-Riley (FR) classes I and II 
which are thought to be the parent populations of respectively FSRQs and BL Lacs. For 
these two classes they derive that the slope of the intrinsic luminosity function is respectively 
2.53±0.06 and 2.65±0.06 which are in good agreement with the value of 2.73±0.38 derived 
here. 

A visual representation of the best fit XLF model (double power law plus MPLE model) 
is shown in the left panel Fig. [10] which reports the volume density of blazar as a function 
luminosity class and redshift. The datapoints are the "deconvolved" BAT observed data, 
that is the number (or density) of blazars which an instrument with optimum sensitivity 
would see. In order to deconvolve the BAT data, we computed for each bin of redshift and 
luminosity, the ratio between the integrals of $(Lx, z) and A(Lx, z) (see Eqs. [2] and [8] for a 
definition of both). This gives a correction factor which allows to deconvolve the BAT data. 
Also note, that given the sparseness of the BAT data, the correction factor is sometimes 
averaged over large bins of redshift and luminosity where the XLF is strongly varying, thus 
it might be somewhat uncertain. Nevertheless, Fig. [TU] highlights that BAT is sampling with 
good accuracy the redshift peak of some of the most luminous objects in the Universe. From 
the same figure it is clear that the density of very luminous blazars (Log Lx > 10^'' erg s~^) 
peaks at large redshift and precisely at z=4.3±0.5. This is much larger than the value of ~ 1.9 



derived (or assumed) for X-ray and Optical surveys (see e.g. lUeda et al.ll2003l : iHasinger et al. 



2OO5I : iBongiorno et al.l 120071 : [Silverman et al.l |2008| ) . The likely reason of this difference will 
be addressed in details in § [61 The right panel of Fig. [10] shows the non-parametric blazar 
XLF built using the 1 /V max method along with the best-fit analytical XLF model (model 
7). It is apparent the good agreement between the two representations. 



4.3. Two populations: FSRQs and BL Lacs 



Previous works (e.g. Wolter et al.lE991 : Rector et ahlboOO : Wolter &: Celotti|[2001 : Beckmann et al. 



2003 



Padovani et al.ll2007D have reported evidence about the different evolutionary behav- 
iors of FSRQs and BL Lacs. The V^/Vmax test reported in § [3] showed that also in our 
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Fig. 10. — Left panel: Number density of blazars (FSRQs and BL Lacs) as a function of 
redshift and luminosity class. The solid lines represent the best-fit XLF model (model 7 in 
Tab. Hj) . The BAT data (points with errors) were "deconvolved" taking into account the 
BAT sensitivity (see § 14.21 for details). Right panel: Luminosity function of the BAT blazars 
built using the 1/Ymax method (datapoints) with superimposed the best fit XLF model 
(model 7 in Tab. H]). 

sample the two classes of objects might evolve differently. In the next sections we test this 
hypothesis. 

4.3.1. FSRQs 

We applied the two best fit models of the previous section (MPLE coupled to a single 
and double power-law local XLF respectively) to the FSRQ class. The best-fit parameters 
are reported in Tab. |H We note that both XLF models produce essentially the same re- 
sult. When the local XLF is modeled as a double power- law model, the faint-end slope 71 
is required to be largely negative (<-50) and the break luminosity coincides with the 
minimum observed luminosity of FSRQs in the BAT sample. Under this conditions, the 
double power-law model reduces to a single power-law distribution with a sharp cut-off at 
Lx < 2 X 10^'' erg s~^. Figures [TT] and [T^ (right panel) show how well the best fit XLF models 
(model 9 and 10 in Tab.Hj) reproduce the observed distributions (in redshift, luminosity and 
source counts). 

Fig. [T^ shows the number, and its volume density, of FSRQs in the Universe for different 
luminosity classes as derived from the best fit. The BAT data were "deconvolved" with the 
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method outlined in § 14. 2[ It is clear that BAT is very effective in constraining the density 
of FSRQs at high luminosity and large redshifts. The same figure shows that the cut-off in 
the evolution is, at least for luminosities larger than 10"^^ erg s~^, well constrained. 




Fig. 11. — Redshift (left) and luminosity (right) distrib ution of the BA T FSRQs. Error bars 
were computed taking into account the Poisson error (jGehrelslll985l ). For both cases, the 
dashed line represents the best fit XLF model (model 10 in Tab. Hj) convolved with the BAT 
sky coverage. 



4.3.2. BL Lacs 

Given the small number of BL Lac objects (12) and the relatively low redshift range 
that they span (0.01 <z< 1.0) we cannot use complex evolutionary models. We thus tried to 
fit a simple PLE model to the data. We obtain an excellent fit which implies mild negative 
evolution (albeit with large errors). Indeed, the best-fit value of the evolutionary parameter 
is -0.79±2.43. As the results reported in Tab. H] show, fixing the evolution parameter at 
zero produces an equally good fit (see model 11 and 12). Fig. [H] shows that the best- 
fit XLF models (model 11) reproduces accurately the observed distributions in redshift and 
luminosity. The best-fit XLF predicts the cumulative source count distribution always within 
1 a (see right panel of Fig. [T2]). We also tried to use a double power-law model for the local 
XLF (see model 13). While this model reproduces the BAT data accurately, most of its 
parameters are poorly constrained. However, we note that the best-fit parameters are in 
good agreement with the values derived for the whole blazars sample (see model 7) and that 
the evolution, although still consistent with zero, became positive. Thus, we believe that 
given the small number of BL Lac objects it is currently impossible to constrain the sign of 
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Fig. 12. — Cumulative log - log S distributions for the BAT FSRQs (left) and for the 
BAT BL Lacs (right). The dashed lines are the predictions from the best fit XLF models 
(model 10 for FSRQs and model 11 for BL Lacs). The different evolution of the two source 
populations is apparent in the different slopes of the two log A^ - log S distributions (steeper 
for FSRQs and flatter for BL Lacs). 




Fig. 13. — All-sky Number (left) and number density (right) of FSRQs as a function of red- 
shift and luminosity class. The solid line is the best-fit XLF model (model 10 in Tab.Hj). The 
BAT data (points with errors) were "deconvolved" taking into account the BAT sensitivity. 



the evolution (i.e. positive or negative evolution). 
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Fig. 14. — Redshift (left) and luminosity (right) distributi on of the BAT BL Lacs. Error 



bars were computed taking into account the Poisson error (jGehrelslll985l ). For both cases, 
the dashed line represents the PLE XLF (single power law model, see model 11 in Tab. Hj) 
convolved with the BAT sky coverage. 

5. Implication for the Cosmic X— ray Background 

Eq. [T7|can be used to estimate the contribution of blazars to the CXB outside the energy 
band of this survey. In this case, Fx{Lx, z) = Fx{Lx, z, E) becomes a function of energy and 
represents the spectral energy distribution (SED) of the blazar source class. To model the 
contribution of FSRQs and BL Lacs we used the best-fit XLF models derived in the previous 
sections. These are model 9 and model 11, in Tab. HI respectively for FSRQs and BL Lacs. 
We also remark that the results presented here do not change if other valid XLF models 
(e.g. model 10 for FSRQs and model 12 and 13 for BL Lacs) are used. In order to compute 
correctly the uncertainties we employed a Monte Carlo simulation. We generated a large 
number (> 200) of luminosity functions starting from randomly sampled best-fit parameters 
drawn from the covariance matrix derived during the fit stage. Moreover, for each randomly 
generated XLF, a random photon index has been drawn from the index distributions of the 
given class (i.e. FSRQ or BL Lac). We then computed the contribution to the CXB for 
each of these luminosity functions and computed the 1 a deviation, around the mean value, 
at given fixed energies. As a first test, we model the SED using a simple power-law model. 
Fig. [T^ shows the contributions of FSRQs and BL Lacs (evolving as different populations) to 
the CXB. It is apparent that while the contribution of BL Lac objects appears negligible in 
this hypothesis (i.e. no BL Lac evolution), the contribution of FSRQs is substantial in hard 
X-rays. From our luminosity function we derive that virtually 100 % of th e CXB for e nergie s 
> 500 keV is produced by FSRQs. We also use the synthesis model of Gilli et al. ( 2007 1^ 
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to take into accoun t the contribution of Seyferts to the CXB. We arbitrarily renormahze 
the iGiUi et aLl (120071 ) model by 1.1. This is justified by the fact that this sy nthesis model is 



tuned to reproduce the CXB as measured by HEAO-1 (IGruber et al 
to underestimate the CXB emission of ~10 % at 30 keV (lAjello et al. 



1999) which is known 



2008b. and references 



therein). It is apparent from Fig. [15] that summing the contribution of blazars to the one of 
Seyferts achieves a good estimate of the intensity of the CXB emission up to the MeV range. 

Modeling the SED with a simple power-law model is a straightforward and robust hy- 
pothesis, but it remains accurate only for extrapolations close to the original 15-55 keV 
band. Indeed, Fig. [15] shows that at 10 MeV the contribution of FSRQs, computed in this 
way, overestimates the diffuse background by an order of magnitude. The uFu spectrum 
of FSRQs exhibits an IC peak which is located somewhere in the MeV-GeV band. While 
detailed modeling of the SEDs of each of the BAT blazars is outside the s cope of this paper 



we note that some of t he BAT FSRQs were analyzed by several authors (IZhang et al. 



2005 



Sambruna et al.ll2006l . 120071 : iTavecchio et al.ll2007l : IWatanabe et al.ll2008l ). In all cases, the 
authors find that the IC peak is located in the MeV band. We thus represent the SED with 
an empirical double power- law model of the type: dN/dE oc [{E/Eb)-^^ + {E/Eb)'^^]'\ 
where Fi and r2 are the photon indices (1.6 and 2.5 respectively) before and after the energy 
break Eb- For the energy break, Eb, we chose a value of IMeV motivated by the observa- 
tions reported above. However, we note that the smooth and large curvature of the model 
we employ, makes it virtually insensitive to the exact value of Eb if this is within 1 order of 
magnitude. Fig. [TH] shows the contribution of FSRQs assuming that their IC peak is located 
in the MeV band. We find that in this case FSRQs account for the entire CXB emission 
up to 10 MeV. While there is basically no difference with respect to the single power-law 
case below 500 keV, the curvature of the IC peak makes the contribution of FSRQs to the 
CXB slightly smaller around 1 MeV. We also note that moving the IC peak beyond 10 MeV 
produces a negligible curvature in the FSRQ integral emission and thus this case is well 
represented by the single power-law model. 

Thus, the two analyses shown here cover well the case in which the IC peak is either 
located at MeV or at GeV energies (double and single power-law model respectively). We 
must therefore conclude that the contribution of FSRQs to the diffuse emission is relevant 
and likely accounts for a substantial fraction (potentially ~100 %) of the CXB around 1 MeV. 
Interpreting the CXB as a strong constraint, we derive that the population of FSRQ sampled 
by BAT must have the IC peak located in the MeV band in order not to overproduce the 
diffuse background at ~10 MeV. iBhattacharya et al.l (120081 ) recently reported for the FSRQs 



■^A web interface to the model of 



Gilli et al. 



(|2007l ) is available at 



http://www.bo.astro.it/~gilli/xrb.litml 
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detected by EGRET a mean photon index of 2.30±0.19. Since FSRQs have a mean photon 
index of 1.6 in BAT, this imphes aheady that the IC peak is located in between the BAT 
and EGRET energy bands. However, as we noted aheady in § [2] only 9 FSRQs are in 
common between the EGRET and the BAT samples and this might imply that the other 
FSRQs detected by BAT have an IC peak at even lower energies. As it w ill discussed i n 
§ 17.21 Fermi-LAT will certainly clarify this scenario. Indeed, very recently, lAbdol (j2009al ). 
discussing the results of the first three months of observations of Fermi-LAT, showed that 
FSRQs are detected by Fermi with a mean photon index of 2.4. Thus, FSRQs have soft 
spectra (photon index > 2.0) in the GeV band while they have hard spectra in the hard 
X-ray band. This confirms that their IC peak is located between the two bands. We note 
th at the shape of t he int egrated emission of FSRQs is similar to the empirically-derived one 
of IComastri et all J2006h . 



We also found that the contribution of BL La cs is very small if they are a non-evolving 



(or mi ldly e volving) popul a tion. In agreement with iGeorgakakis et al.l (120041 ) , iGalbiati et al. 



(120051 ) and iGiommi et al.l (120061 ) we find that the contribution of blazars (FSRQs and BL 
Lacs) to the 2-10 keV CXB is ~10%. 



6. The high-redshift non-thermal Universe 



The fact that the shape of the blazar luminosity function and its evolution are in agree- 
ment with those of X-ray selected AGNs sug gests the presence of a link between accre- 
tion and jet activity (e.g. iMerloni et al.l l2003l ). In other words, it seems that the most 
luminous AGNs (which in turn are the most luminous QSOs) harbor the most power- 
ful blazars. This scenario takes place mostly in the very high-redshift Universe where, 
thanks to the abundance of dust and gas, efficient accretion led to the build-up of mas- 
sive QSOs. However, their space density quickly decreases, with cosmic time, leaving the 
room for the bulk of low-luminosity QSOs. This "anti-hierarchical" scer iario, where larger 



structures conie firs t, was also named "cosmological downsizing" (e.g. iCowie et al.l Il999 
Hasinger et al.l l2005l . and references therein) and constitutes a unique phenomenon which 
is not pred icted in most of the serni-ana l ytic models based on Cold Dark Matter structure 
formation (IKauffmann fc Haehneltl l2000l : IWyithe fc Loebl 120031 ) . The late evolution of low- 
luminosity AGNs coincid es well with the peak of the star formation in the Universe (e.g. 
Hopkins fc BeacomI |2006| ) highlighting once more the interconnection between the host and 
its nucleus. A mechanism of accretion with different efficiencies, as a fun ction of cosmi c 
time, has been invoked to explain the anti- hierarchical growth of AGNs (e.g. iMerlonil 120041 ). 



However, the fate of the very first and luminous quasars (i.e. the apparent disappearance of 
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Fig. 15. — Spectrum of the CXB and contribution of the FSRQs (blue region). The 
dat a points are different measurements of the diffuse background as indicated in the la- 



bel flFukada et al. 



2000; 



1975 



Revnivtsev et al. 



Gendreau et al.lll995l: IWatanabe et al.lll997l : IWeidenspointner et al. 



2OO3I : kiello et al.ll2008bh . The dashed line is the total contribution 



of Seyfert-like AGNs computed with the model of iGilli et al.l (120071 ) arbitrarily multiplied 
by 1.1 to fit the CXB emission at 30keV. The solid line is the sum of the Seyfert-like and 
FSRQs. The spectrum of FSRQs has been modeled as a power-law with a mean photon 
index of 1.6. The blue region represents the range of values obtained from the Monte Carlo 
realizations of best-fit parameter ranges. The magenta solid line represents the contribution 
of BL Lac objects whose uncertainty is not plotted for clarity, but is, due to the low number 
of objects, > 30 % at any energy. 



quasar activity in massive galaxies at late times) remains stil l unknown and there are doubts 
whether these objects can form at all in a ACDM Universe (jSpringel et al.ll2005l ). 
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Fig. 16. — Contribution of FSRQs (blue region) to the CXB. The data are the same as in 
Fig. [151 but in this case the SED of the FSRQs has been modeled with a double power-law 
function. The IC peak is located in the ~MeV region. The contribution of BL Lacs is the 
same as in Fig. [15] and is not drawn here for clarity. The blue region represents the range of 
values obtained from the Monte Carlo realizations of best-fit parameter ranges. 



In Fig.[T7]we compare the shape of the evolution o f the most luminous BAT blazars with: 
1) th e evolution of luminous X-ray selected AGNs (IHasinger et all I2OO5I: [Silverman et al.l 



20081 ). 2) the star formation history (SF R) of the Universe ([Hopkins fc Beacom 



200 



and the evolution of UV and IR galaxies ( Bouwens et al. 20081 : ISandersI 2004 . respectively) 
The most direct comparison is clearly with AGNs. We note that the shape of the evolution 
is very similar, but that the cut-off in the AGNs growth is at lower redshifts (z<2) with 
respect to the cut-off of BAT blazars. This means that the peak of the evolution of the BAT 
objects is at much earlier times in the history of the Universe than the peak of the most 
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luminous AGNs detected in the deepest X-ray surveys. The star formation history of the 
Universe shows a similar trend and a peak around z=~2 as for normal AGNs and this has 
been interpreted as the evidence of the strong link between AGNs and the star formation 
in its host galaxy. Assuming that the activity of the BAT blazars is powered by accretion 
onto super-massive black hole, this implies that the doppler-boosting allows BAT to detect 
a class of objects (in luminosity) which escaped even the deepest X-ray surveys. 




Fig. 17. — Density of the most luminous BAT blazars (47.3<Log Lx <48.3, blue datapoints 
and line) as a function of redshift compar ed to: densities of AGNs from X-ray surveys 
(IHasinger et al.ll2005l : [Silverman et al.ll2008h. star formation rate (jHopkins fc Beacomll2006l ) 
and densities of UV and IR galaxies ( Bouwens et al. 2008 : SandersI 12004 . respectively). The 
different curves were rescaled arbitrarily to match the evolution of the BAT blazars. Thus, 
the comparison involves only the shape. 



The only objects which show a peak in the evolution at redshift >2 are bright star- 
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forming galaxies det ected up to very high redshifts in the GOODS fields (IWall et al.l 12008 
Bouwens et al.ll2008l . and references therein). The rapid brightening of galaxies within the 
first two billion year s is ti ghtly connected to the assembly of large dark matter haloes (e.g. 
Wang &: KauffmannI l2008l ). An intriguing idea is that black holes are formed and fueled , 



and AGN activity is triggered during maj or galaxy mergers (IKauffmann fc Haehneltl 12000 



Wyithe &: Loebl l2003l : ICroton et al.l 120061 ) . Although this constitutes a plausible fueling 
mechanism, the interr nediate- lumin osity AGNs are harbored by galaxies which display little 
or no merger events (IHasanl 120071 ). However, to form a super- massive black hole, a more 
violent process such as a major merger event may be required to funnel a large amount of gas 
into the central region of the galaxy. It seems that the most massive galaxies undergo a major 
merger event at earlier times, and within the fir st 2-4 billion years, than less massive galaxies 
( iWang fc Kauffmannll2008l : IStewart et al.ll2008l ). If the bulk of the black hole mass is formed 
in this way, then it would explain the lack of growth of powerful AGNs at present times. 
The abundance of luminous blazars at large redshift fits well in this scenario since blazars 
are found in gian t elliptical galaxies which in turn are supposed to have undergone a ma- 
jor rn erger event (IToomre fc Toomrd 1 19721 : iNegroponte fc White! Il983l : IWang fc Kauffmann 
20081 ). Thus, we argue that the luminous blazars can be used as a tracer of massive galaxies 
and merging activity in the very early Universe. 

Understanding the formation of massive galaxies is an important astrophysical issue 
because as mu ch as 50 % of th e stellar mass in the local Universe appears to be in early- 
type systems ( iBell et al.l l2003l ). The most massive galaxies (i.e. > 10^^'^ Mq) appear 
to be already i n place at z~2 suggesting that they formed in the very early Universe (see 
Conselicd l2008l . for a review). Unfortunately present surveys, both in the optical and in 
X-rays, are not sensitive enough to make a statistical census beyond redshift 2. Here, BAT 
and more in general blazar surveys, can play an important role. In Fig. [18] we compare 
the redshift evolution of the luminous BAT blazars with the p rediction of the evolut ion of 
massive elliptical galaxies as determined, using simulations, by iDe Lucia et al.l (120061 ). The 
similarity between the two curves is apparent and reinforces our idea that blazars can be 
used to study the formation of massive systems in the early Universe. 

The number density of massive ellipticals hosting active blazars is larger than the 
one reported in Fig. [TS] since we consider only those blazars pointing at us. Since there 
should be no difference between the population of blazars pointing at us and those point- 
ing in all the other directions we might try to estimate their total number. The number 
of misali igned blazars, can be estima ted as: Ntot ~ 2r^, where F is the bulk Lorentz fac- 
tor (e.g. 



Scheuer fc Readhead 



1979[). Assuming a standard value of F = 10 (20) (see e.g. 



Urry &: Padovani IQQsl : Sambruna et al. 2007) leads to a space density of ellipticals, hosting 
active blazars, which is 200 (800) times larger than that one shown in Fig. [18] and thus of 
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Fig. 18. — Density of the most luminous BAT blazars (47.3<Log Lx <48.3, blue data- 
points and line) comp ared to the pre diction of the evolutio n of massive elliptical galaxies of 
De Lucia et al.l (120061 ). The curve of iDe Lucia et al.l (120061 ) was rescaled to match the BAT 
luminosity function. In order to consider also the contribution of misaligned jets, the space 
density of blazars (y-axis) needs to be multiphed by a factor in the 200-800 range (see § [6] 
for details). 



the order of 2.0-8.0x 10"^ Mpc'^ at a redshift of ~4. 

Moreover, the excellent agreement between the evolution of massive galaxies and blazars 
suggests that there should be no intrinsic difference, at least in the early Universe, between 
ellipticals and ellipticals hosting active blazars. This means that the fraction of ellipticals 
which host active blazars is not changing dramatically as a function of cosmic time. 



Recently, 



Sikora et al, 



(120071 ) showed that on a radio-loudness/Eddington-ratio dia- 
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gra m ellipticals and spirals form two distinct and well-separated parallel sequences (see 
also IWilson fc ColbertI Il995l ) . They argue that this different behaviour might be given by 
the spin. Indeed, within the hierarchical cosmological framework, the main difference in 
the evolution of giants ellipticals and spirals is that the first ones underwent at least one 
major merger event in the past. These events can produce a maximally rotating b lack hole 
by coalescence of the two black holes (lEscala et al.l 12004 l2005l: iDotti et al.l 120071) and by 



trigg e ring large-gas accretion events that spin up the hole (IBarnes fc Hernquistlll996l : lEscala 



20071 : iBogdanovic et al.l 120071 ). It is import ant to note that und e r the assumption that jets 



are powered by rotating black holes via the iBlandford fc Znajeki (119771 ) mechanism, the effi- 
ciency of the jet production is determined by the black hole spin. We believe that gathering 
a large sample of blazars which spans adequately luminosity and redshift will allow in detail 
to understand the role of spin and accretion in triggering jet activity. In this respect, BAT 
and, in particular, Fermi will play an important role. 



7. Summary and Conclusions 



7.1. The blazar XLF 



We have used a complete sample of blazars detected by BAT to derive the first luminosity 
function of blazars in the 15-55 keV band. We have shown using several methods that BAT 
blazars are evolving strongly while the Seyfert-like AGNs detected by BAT are not. The 
evolution of the blazars is implicit in the fact that BAT detects 10 objects at redshifts larger 
than 2.0. The local luminosity function (e.g. ^{Lx,z = 0)) is compatible with a double 
power-law model where the faint-end sl ope is required t o be f iat mainly by the absence of low- 
luminosity BAT blazars. According to lUrry fc Shaferl (jl984j ) the flattening of the local XLF 
at low luminosities might be produced by beaming which boosts intrinsically low-luminosity 
sources to high luminosities. The best-flt XLF models imply an evolution in luminosity 
which is epoch dependent (e.g. the evolutio n parameter changes with redshi ft) as was found 



for other samples of blazars or FSRQs (e.g. IWallll2008l : iPadovani et al.ll2007l ). In all cases, it 



appears that blazars are evolving strongly up to a redshift cut-off which is, at least for the 
most luminous objects, well constrained by our data. 

We flnd that the strong evolution of blazars is driven by the evolution of FSRQs which 
are the only objects detected by BAT at large redshifts. Our best-flt XLF, to the sample of 
BL Lacs, shows tha t BL Lacs ha v e a ne g hgible evolution . This is found to be in agreement 



with the results of iRector et al.l (120001 ) , iPadovani et al.l (120071 ) and of iBhattacharya et al. 



( I2OO8I ). Given the small number of objects in o ur sample we cannot rule out, nor conflrrn . 



the claims of negative evolution of BL Lacs (e.g. IWolter et al.lll99ll : iBeckmann et al.ll2003l ). 
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These last two samples reach a flux which is lower than the current BAT sensitivity, thus 
with a few more years of exposure BAT might be able to test this negative-evolution scenario. 
Finally, we also remark that given the small number of objects (12), the evolution of BL 
Lacs is marginally consistent (at 1.5 o") with the evolution of the FSRQ class. 

The log A^- logS" distribution of blazars is steeper than the one of Seyfert-like AGNs. 
Its slope of 1.9 is larger than the Euclidean value of 1.5 which characterizes the Seyfert 
galaxies (see Tab. [3] for details). We, thus, expect that the fraction of blazars will steadily 
increase among the total AGN population detected by BAT. We expect that in a relative 
short timescale (e.g. a couple of years) and, depending also on the systematic errors of 
the BAT survey, the blazar sample might contain more than 60 objects. This will be very 
important as it will allow to improve the results and the prediction of the XLF models for 
both FSRQ and BL Lac objects. As we have shown, blazars are extremely rare objects 
which can be detected in hard X-rays only through large-area sky surveys. Our best-fit XLF 
model predicts a flattening i n the log A^-lo g S at fluxes lower than 5x10^^^ erg cm^^ s^^. 



Thus a mission like EXIST (lGrindlayll2005l ) would detect FSRQs with a surface density of 



~0.5deg~^ at fluxes of 5xl0~^^erg cm~^ s~^. If confirmed that most of the BAT FRSQs 
are MeV blazars (see § IHl and UTI^i . then an optimum band to s elect and study them would 
be the MeV band. A mission like GRIPS ( IGreiner et al.ll2008l ) would gather a fairly large 
sample (500-1000) of blazars. 

The redshift distribution of the BAT blazars (e.g. Fig. [S]) shows a peak at low redshift 
and a flat tail extending up to zp^ 4. This distribution differs from the redshift distributions 



of radio-selected blazars which display a peak at zf»1.0-1.5 (e.g. iDunlop fc Peacockl Il990 



Wall et al.l |2005| ) . The reasons for this difference lie in the different selection effects and 
sensitivity of these surveys. Since BAT is sensitive only to bright X-ray fluxes most of the 
low-luminosity low-redshift sources are currently undetected. Additionally, the discrepancy 
in the redshift distributions between radio and hard X-rays might also be due to the different 
shapes of the evolution and of the local luminosity functions in these bands. This would not 
be surprising in view of the fact that radio and X-rays probe different scales in these systems. 
A larger dataset of hard X-ray-selected blazars will allow us to test these hypothesis. 



From our XLF we derive that the density of FSRQs at fluxes of 10 er g cm ^ s 



5.2^^^ deg ^ which is compatible wit h previous estimates in other X-ray bands (IWolter fc Celotti 



IS 



200ll : iGiommi fc Colafrancescdl2006l ). Within the large uncertainties of our estimate, we de- 
rive that at faint fluxes the density of FSRQs is not negligible when compared to the total 
AGN population. As an example, the densi ty of all AGNs in the XMM-COSMQS field is 
24.0=b3deg~^ for equivalent fluxes as above (ICappelluti et al.l 120071 ). This means that deep 
X-ray surveys necessarily contain a ~10 % fraction of blazars. This seems in agreement with 
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the finding of Idella Ceca et al.l (119941 ) wlio report a fraction of radio-loud objects (among X- 



ray selected AGNs) of ~1 % at fluxes of 10 erg cm ^ s ^. Simi l ar fra ctions of radio-loud 



AGNs were also found by iHooper et al.l (Il996l ) and IZickgraf et al.l (120031 ) . 



The main point of concern is however the selection of very absorbed, Compton-thick, 



Fiore et al. 


2008; 


Brusa et al. 


2008) 



The current 'paradigm' is that exceptionally hard X-ray spectra (photon indices of 1.0-1.5) 
are likely produced by strong absorption. Our analysis shows that FSRQs have intrinsically 
hard X-ray spectra with photon indices so metimes lower than 1.6 (see Fig. [21 but also 



Tavecchio et al.l 120071 : IWatanabe et al.ll2008l ). We believe that, if radio properties are not 
properly taken into account, selection of absorbed sources based solely on hardness ratios 
will produce a sample which can be contaminated by a substantial fraction of FSRQs. 



7.2. The Cosmic X-ray Background 



The origin of the MeV background has been a long-standing issue in astrophysics. Sev- 
eral astrophy sical processes have be en put forward to explain it. Among ther n, dark matter 
annih ilation ( Ahn fc Komatsulboosl ) . nuclear decays fro m Type la superno vae ( Clayton &: Ward 
1975 ) and non-thermal emission from Seyfert galaxies ( Inoue et al.lboosi ) were the most im- 
portant ones. We used our best-fit XLF model to make a prediction of the integrated emission 
due to FSRQs and derived that FSRQs account for most of the diffuse background emission 
for energies >500keV. Moreover, assuming that most of the FSRQs have an IC peak in 
the MeV band, as some of the BAT blazars (e.g. IZhang et al.ll2005l : ISambruna et al.ll2007l : 



Tavecchio et al.ll2007l : 



Watanabe et a. 



20081 ). we showed that the sum of the contribution of 



emission-line AGNs ( IGilli et al.l 120071 ) and blazars reproduces well the CXB emission from 
1 keV to 10 MeV. Our prediction of the contribution of blazars to the CXB is well in agree- 
ment, in the 2-10 keV band, with the fin dings of several authors (e.g. 
Comastri et al.ll2006l : iGiommi et al.ll2006l ). 



Galbiati et al. 2005 



Recently llnoue et al.l (120081 ) proposed that a population of non-thermal electrons present 
in the hot AGN coronae can account for a substantial part of the MeV background. Our 
finding shows that the non-thermal contribution from AGN coronae should be small as 
most of the diffuse ba ckground emission is accounted for by blazars. In a more recent work, 
Inoue fc Tot anil (120081 ) derived the luminosity function of EGRET blazars taking into account 
the blazar sequence; this is then used to compute the contribution of blazars to the diffuse 
background. From their best fit model, it arises that blazars contribute negligible emission 
around lOkeV. This is in conflict with the main finding of this paper that blazars contribute 
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~10-20 9^ of the CXB emission in the 15-55 keV band. As a matter of fact, ~17% of all 
BAT AGNs are blazars and thus their contribution to the CXB must be of the same order. 

The scenario which we derive from our data can be easily tested by the Fermz-Large Area 
Telescope (LAT). Indeed, we showed that in order not to overproduce the MeV background, 
most FSRQs are required to 'peak' at MeV energies for a large fraction of their time. Thus, 
the detection by LAT of soft FSRQs (e.g. photon indices of 2.2-2.5) would constitute a 
final evidence that the IC peak should be between the BAT and the LAT energy bands. 
On the other hand, the detection of hard FSRQs (indices of 1.4-1. 8) would inva lidate our 



prediction. The first Fermi-L AT results convalidate our results (lAbdol l2009al ). Indeed, 
FSRQs are detected by Fermi with a mean photon index of 2.4 (and a tail extending up to 
3.0) confirming that the IC peak must be located somewhere between the keV and the GeV 
band. 



7.3. Tracing the star formation history of massive ellipticals at high redshift 



The similar evolution of radio-quiet AGNs and star formation history of normal galaxies 
has been interpreted as the evidence of the co-evolution of AGNs and their hosts (e.g. 
Madau et all 1 19961 : iHasinger et aDbnn,'^! ). 



The main, serendipitous, finding of our analysis is that the evolution of the BAT blazars 
shows a redshift cut-off which is larger than previously found for other, mostly radio-quiet, 
AGN samples. This is found to be Zc =4.3±0.5 for blazars of typical luminosities exceeding 
10^^ erg s^^ . The large redshift cut-off shows that the most luminous blazars formed very 
early in the Universe and then their number density quickly decreased. To our knowledge, no 
other source class displays a simila r extreme evolution. X-ray surveys show that the redshift 



cut-off increases with lumi nosity (lUeda et al.l l2003l : iHasinger et al.l l2005l : iLa Franca et al. 



20051 : [Silverman et al.ll2008l ) and thus we believe that doppler-boosting, due to the relativistic 
beaming, allows BAT to detect rare objects which escaped even the deepest surveys. 

We compared the blazars lumin osity function and th e prediction of the star formation 
history of massive elliptical galaxies (jPe Lucia et al.ll2006l ) and found good agreement. How- 
ever, this agreement is not entirely surprising if one realizes that blazars are normally foun d 
in giant elliptical galaxies (lUrry fc Padovanilll995l : iFalomo et al.ll2000l : lO'Dowd et al.ll2002l ). 
This represents another evidence that AGNs (jet activity in this case) and their hosts co- 



^The BAT blazar XLF's main uncertainty is given by the small number of objects. Moreover, the CXB 
fraction increases if FSRQs and BL Lacs are treated as a single population. 
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evolve through the history of the Universe. However, tracing the evolution of giant galaxies 
is currently at the limit, or beyond, of the present-generation instruments and thus the use of 
blazars might represent the only approach to understand the formation of the most massive 
galaxies in the early Universe. 

Elliptical galaxies are thou ght to be the only objects which undergo one major merger 



[e.g. iWang fc Kauffmanru |2008| ) and in particular this seems to happen in the first bil- 
lion years of the Universe. As a natural consequen ce, it is believed tha t merging activ- 
ity would produce a rapidly spinning black hole (e.g. IVolonteri et al.l 120071 ) which on theo- 
retical grounds, is required to explain the produc ti on of a coUimated, relativistic, outflow 
Jsiandford fc Znaiekl [l977l : Isiandford et~al] Il990h . Isikora et all J2007h found out that on 
a radio-loudness/Eddington-ratio diagram elliptical and disk/spiral galaxies form different 
sequences and invoke the spin as the black hole parameter which might explain this different 
behavior. Larger blazar samples, better understanding of the evolution of massive systems, 
and direct black hole spin measurements will help in clarifying the jet-spin- merger scenario. 
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Table 3. Y/Vmax and log N- log S tests. 



Sample 


<V/Vmax> 


(3^ 


Ab 


# Objects 


Scyfcrts 


0.509±0.021 


1.496±0.073 


6.70±0.48 


199 


BLAZARs 


0.666±0.045 


1.932±0.206 


1.27±0.20 


38 


FSRQs 


0.728±0.056 


2.077±0.269 


0.83±0.16 


26 


BL Lacs 


0.576±0.083 


1.694±0.316 


0.38±0.10 


12 



^Best-fit exponent of the log A^-log S distribution (e.g. A'"(> S) = 
A S-f^). 

^Surface density of objects above 10~^^ erg cm~^ s~^ in units of 
10-3 deg-2. 



Table 4. Parameters of fitted Luminosity Functions. Parameters without an error estimate were kept fixed during 
the fitting stage. Kolmogorov-Smirnov (KS) values are the probabilities of the model and the data to be drawn from 

the same parent population. 



Octill jJiU 


44- Ohifrf^ 

■ff- WDJCljlO 


MoHpI 

IVlOtiCl -ff- 






71 


■^9 
7Z 




]j 




7 






CXB % 


Seyferts 


199 


1 


PLE+2pow 


0.909+0.064'= 


0.80+0.08 


2.67+0.20 


0.61+0.18 









0.99 


0.94 


21.6% 


Seyferts 


199 


2 


PLE+2pow 


0.778+0.055'' 


0.84+0.08 


3.01+0.30 


0.61+0.14 


2.62+1.18 






0.99 


0.93 


55.2% 


BLAZARc 


38 


3 


PDE+lpow 


0.757+0.135 




2.67+0.13 


1.0 


4.00+0.77 






0.43 


0.22 


22.0% 


BLAZAR"= 


38 


4 


MPDE+lpow 


0.732+0.117 




3.08±0.20 


1.0 


8.95+1.90 




-0.69±0.27 


0.34 


0.40 


318.0% 


BLAZAR': 


38 


5 


MPLE+lpow 


0.804+0.131 




3.13+0.21 


1.0 


2.96+0.47 




-0.23+0.08 


0.43 


0.41 


18.4 % 


BLAZAR 


38 


6 


PDE+lpow 


0.255+0.041 




2.26+0.07 


1.0 


2.05+0.57 






0.002 


0.00 


6.9% 


BLAZAR 


38 


7 


MPLE+2pow 


1.379+0.224 


-0.87+1.31 


2.73+0.38 


1.81+0.77 


3.45+0.44 




-0.25+0.07 


0.86 


0.88 


~20.0% 


BLAZAR 


38 


8 


MPDE+2pow 


0.948+0.152 


-0.83+1.43 


2.54+0.21 


1.95+0.93 


11.62+1.40 




-0.85+0.18 


0.46 


0.85 


640.0 % 


FSRQd 


26 


9 


MPLE+lpow 


0.533+0.104 




3.45+0.20 


1.0 


3.72+0.50 




-0.32+0.08 


0.86 


0.87 


9.0% 


FSRQ" 


26 


10 


MPLE+2pow 


0.175+0.034 


< -50.0 


2.49+0.37 


2.42+0.19 


3.67+0.48 




-0.30+0.08 


0.85 


0.89 


8.3% 


BLLac'^ 


12 


11 


PLE+lpow 


0.830+0.240 




2.61+0.37 


1.0 


-0.79+2.43 






0.55 


0.33 


0.3% 


BLLac": 


12 


12 


PLE+lpow 


0.784+0.226 




2.73+0.17 


1.0 









0.55 


0.33 


(!).3% 


BLLac 


12 


13 


PLE+2pow 


1.506+0.435 


-0.89+3.7 


2.51+1.61 


1.78+2.71 


1.54+2.73 






0.71 


0.93 


1^2% 



I 



''In unit of 10 ^ Mpc ^ erg ^ s unless otherwise stated. 
''In unit of 10-^ Mpc-^. 

LMin has been set to 6 x 10**^ erg~' s. 
'^I^Min has been set to 2x10'*'' erg"-' s. 

''Given the very small value of faint-end slope 71, this model is equivalent to a model where the local XLF is parametrized as a single power law with a sharp cut-off at the lowest 
observed FSRQ luminosity of ~ 2 x lO'*" erg s"'. 
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A. The Seyfert Sample 



In this section we report the 199 Seyfert objects which are used as a control sample 
plus the 6 radio galaxies detected by BAT. The sample, shown in Tab. [5], is reported here 
as a reference for the reader to demonstrate that the classifications reported in Tab. [2] are 
accurate. We note that 152 out of the 205 sources reported in this table are also detected 
in the BAT 22 month survey of iTueller et al.l (120091 ). The main differences among the two 
analyses are: 



th e different ene r gy ba nd used (15-55 keV band versus the 15-195 keV band adopted 
bv lTueller et~aD fl2009h ). 



the di fferent exposure used (36 months versus the 22 months used by iTueller et al. 

jiogl)). 



the slightly differe nt data filtering and screening techniques (see lAjello et al.l l2008a 



Tueller et al.ll2009l . for details) 



Despite these di fferences, ~75 % of the sources detected in this analysis are also contained 
in the sample of ITueller et al.l (120091 ). 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 
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Table 5. Sample of Seyferts and Radio Galaxies 



SWIFT NAME 


R.A. 
(J2000) 


Decl. 
(J2000) 


Flux 
(10-" cgs) 


S/N 


ID'' 


Typeb 


Redshift 


In BAT 22 months^ 


J0006.4+2009 


1.600 


20.152 


1.16+0.20 


5.8 


Mrk 335 


Syl 


0.0254 


y 


J0038.6+2336 


9.650 


23.600 


1.10+0.21 


5.3 


Mrk 344 


Sy 


0.0240 




J0042. 7-2332 


10.680 


-23.548 


2.44+0.21 


11.7 


NGC 235A 


Sy2 


0.0222 


y 


J0048.7-I-3157 


12.188 


31.962 


7.71+0.20 


37.8 


Mrk 348 


Sy2 


0.0150 


y 


J0051.9-I-1726 


12.998 


17.447 


1.81+0.21 


8.6 


QSO B0049+171 


Syl 


0.0642 


y 


J0059.9+3149 


14.997 


31.831 


1.66+0.21 


8.0 


SWIFT J0059.4+3150 


Syl 


0.0149 


y 


JOIOI. 0-4748 


15.274 


-47.800 


0.97+0.18 


5.6 


2MASX J01003469-478303 


GALAXY 


0.0753 


y 


J0108. 8+1321 


17.201 


13.351 


1.78+0.22 


8.2 


4C 13.07 


Sy2 


0.0596 


y 


JOlll. 4-3805 


17.867 


-38.086 


1.52+0.18 


8.3 


NGC 424 


Sy2 


0.0116 




J0113. 8-1450 


18.453 


-14.850 


1.24+0.21 


5.8 


Mrk 1152 


Syl 


0.0522 


y 


J0114. 3-5524 


18.600 


-55.400 


0.92+0.17 


5.3 


SWIFT J0114.4-5522 


Sy2 


0.0121 




J0123. 8-5847 


20.952 


-58.785 


2.65+0.17 


15.3 


Fairall 9 


Syl 


0.0470 


y 


J0123. 8-3501 


20.974 


-35.067 


2.72+0.18 


14.7 


NGC 526A 


Syl.5 


0.0191 


y 


J0127.9-1850 


22.000 


-18.847 


1.27+0.20 


6.2 


MCG-03-04-072 


Syl 


0.0430 




J0134.0-3629 


23.506 


-36.486 


2.36+0.18 


13.0 


NGC 612 


GALAXY 


0.0298 


y 


J0138.6-4000 


24.674 


-40.008 


3.17+0.18 


18.0 


ESQ 297-018 


Sy2 


0.0252 


y 


J0142.6-I-0118 


25.652 


1.300 


1.28+0.22 


5.7 


[W2003c] J014214.0+011615 


Syl 


0.0500 




J0152.9-0326 


28.250 


-3.448 


1.47+0.22 


6.6 


ICR J01528-0326 


Sy2 


0.0172 


y 


J0201. 2-0649 


30.320 


-6.821 


4.17+0.22 


19.3 


NGC 788 


Sy2 


0.0136 


y 


J0206. 5-0016 


31.631 


-0.270 


1.53+0.22 


6.9 


MRK 1018 


Syl.5 


0.0424 


y 


J0215. 0-0044 


33.751 


-0.749 


1.30+0.22 


5.9 


Mrk 590 


Syl.2 


0.0265 




J0226.0-6315 


36.500 


-63.250 


0.91+0.18 


5.2 


FAIRALL 0926 


Syl 


0.0580 




J0226. 8-2819 


36.703 


-28.324 


1.14+0.18 


6.4 


2MASX J02262568-2820588 


Syl 


0.0600 




J0228.4-I-3118 


37.120 


31.316 


4.38+0.23 


19.4 


NGC 931 


Syl.5 


0.0166 


y 


J0232.0-3639 


38.020 


-36.662 


1.09+0.17 


6.4 


IC 1816 


Sy2 


0.0169 


y 


J0234.4-I-3229 


38.612 


32.489 


1.60+0.23 


7.1 


NGC 973 


Sy2 


0.0167 


y 


J0234.8-0847 


38.702 


-8.794 


2.13+0.21 


10.2 


NGC 985 


Syl 


0.0430 


y 


J0235. 6-2935 


38.900 


-29.600 


0.99+0.18 


5.6 


ESQ 0416-G0002 


Syl.9 


0.0592 


y 


J0238. 5-5213 


39.647 


-52.220 


1.31+0.17 


7.6 


ESQ 198-024 


Syl 


0.0452 


y 


J0239. 0-4043 


39.767 


-40.732 


0.97+0.17 


5.8 


2MASX J02384897-4038377 


Syl 


0.0610 




J0241.5-0813 


40.381 


-8.220 


1.34+0.21 


6.4 


NGC 1052 


Sy2 


0.0050 


y 


J0242.9-0000 


40.732 


-0.012 


2.00+0.22 


8.9 


NGC 1068 


Sy2 


0.0038 


y 


J0249.3-I-2627 


42.349 


26.451 


1.25+0.23 


5.5 


IRAS 02461+2618 


GALAXY 


0.0580 


y 


J0252. 8-0830 


43.200 


-8.500 


1.06+0.21 


5.0 


MCG-02-08-014 


Sy2 


0.0168 


y 


J0255.4-0010 


43.873 


-0.170 


4.48+0.22 


20.1 


NGC 1142 


Sy2 


0.0288 


y 


J0256.4^3212 


44.117 


-32.208 


1.31+0.17 


7.7 


ESQ 417-6 


Sy2 


0.0164 


y 


J0311.6-2045 


47.919 


-20.760 


1.27+0.18 


6.9 


2MASX J03111883-2046184 


Syl 


0.0660 


y 


J0325. 1-1-3409 


51.296 


34.152 


1.61+0.24 


6.8 


2MASX J03244119+3410459 


Syl 


0.0629 


y 


J0333. 5+3716 


53.397 


37.278 


1.63+0.24 


6.8 


ICR J03334+3718 


Syl 


0.0574 




J0333. 7-3608 


53.433 


-36.141 


3.12+0.17 


18.9 


NGC 1365 


Syl.8 


0.0055 


y 


J0342. 2-2114 


55.554 


-21.244 


2.15+0.18 


11.8 


SWIFT J0342.0-2115 


Syl 


0.0145 


y 


J0347.3-3029 


56.850 


-30.500 


0.89+0.17 


5.3 


RBS 0741 


Syl 


0.0950 




J0350.7-5022 


57.679 


-50.377 


1.29+0.17 


7.5 


SWIFT J0350. 1-5019 


GALAXY 


0.0365 


y 


J0357.0-4039 


59.268 


-40.666 


0.89+0.17 


5.4 


2MASX J03565655-4041453 


GALAXY 


0.0747 




J0402.5-1804 


60.639 


-18.077 


1.35+0.19 


7.0 


ESQ 549- G049 


Sy2 


0.0262 


y 
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Table 5 — Continued 



SWIFT NAME 


R.A. 

(J2000) 


Decl. 
(J2000) 


Flux 
(10-" cgs) 


S/N 




Typeb 


Redshift 


In BAT 22 months'^ 


J0407. 5+0342 


61.883 


3.717 


1.90+0.25 


7.6 


3C 105 


Sy2 


0.0890 


y 


J0415. 2-0753 


63.800 


-7.900 


1.31+0.23 


5.6 


LEDA 14727 


Syl 


0.0379 


y 


J0426. 4-5712 


66.603 


-57.201 


1.40+0.17 


8.2 


IH 0419-577 


Syl 


0.1040 


y 


J0433. 4-1-0521 


68.355 


5.365 


5.21+0.26 


19.8 


3C-120 


Syl 


0.0330 


y 


J0438.5-1049 


69.633 


-10.830 


1.48+0.23 


6.4 


MCG-02-12-050 


Syl 


0.0360 


y 


Tf\ A A A ^ nolo 

J0444. 7-2812 


71.199 


-28.200 


1.07+0.18 


5.9 


0"R 4" A Cl ■\7" T f\ A A ^ 000000>< 

2MAt3X J 04450628-2820284 


Sy2 


0.1470 




J0451. 8-5807 


72.966 


-58.133 


0.88+0.17 


5.2 


RBS 0594 


Syl 


0.0900 


y 


J0453. 5+0403 


73.380 


4.060 


2.11+0.28 


7.6 


CGCG 420-015 


Sy2 


0.0296 


y 


J0455. 3-7528 


73.841 


-75.477 


1.27+0.18 


6.9 


ESO 33-2 


Sy2 


0.0184 


y 


J0505. 9-2351 


76.497 


-23.854 


2.78+0.20 


13.9 


XSS J05054-2348 


Sy2 


0.0350 


y 


J051d. 2-0009 


79.071 


-0.161 


4.11+0.28 


14.5 


/^O/^ T*or"i o ooo 

(4bO-B0513-002 


Syl 


0.0327 


y 


J 05 19.7-3240 


79.930 


-32.676 


2.38+0.19 


12.9 


sWll' 1 J0519. 5-3140 


Sy2 


0.0350 


y 


J0519.8-4546 


79.963 


-45.774 


2.49+0.17 


14.9 


Pictor-A 


Syl 


0.0351 


y 


J0524. 2-1212 


81.050 


-12.200 


1.42+0.25 


5.6 


T T\ A 1 ^O O O 

LEDA 17233 


Syl 


f\ f\ A OO 

0.0490 


y 


J0552. 3-0727 


88.090 


-7.457 


14.75+0.29 


51.7 


NGC 2110 


Sy2 


0.0078 


y 


J0552. 3+5929 


88.100 


59.500 


"1 "1/1 1 ri "1 

1.14+0.21 


5.3 


IHAb 05480+597 


Syl 


0.0585 




J0558. 1-3820 


89.549 


-38.347 


2.12+0.18 


11.6 


EXO 055620-3820.2 


Syl 


0.0340 


y 


J0559. 9-5026 


89.980 


-50.441 


0.94+0.17 


5.7 


PKS 0558-504 


RG 


0.1370 




J0o02. 9-8633 


90.749 


-86.555 


1.82+0.22 


8.4 


bWll" 1 J0601. 9-8636 


Sy2 


0.0064 


y 


J0d03. 1+6523 


90.799 


65.399 


1 oo 1 n on 

1.38+0.20 


6.8 


UGO 3386 


ATA VV 


0.0154 




J0615. 8+7101 


93.967 


71.021 


6.08+0.20 


30.8 


Mrk 3 


Sy2 


0.0135 


y 


J0d23. 9-3214 


95.994 


-32.248 


1.53+0.20 


7.5 


EbO 426-G 002 


/~1 ATA "VA/" 

GALAXY 


0.0224 


y 


J0624. 1-6059 


96.028 


-60.998 


1.25+0.17 


7.4 


SWIFT J2141. 0+1603 


Sy2 


0.0410 




J0640. 7-4324 


100.200 


-43.400 


r\ oo 1 o 1 o 

0.92+0.18 


5.2 


OA T A O "V T0/^><00/^00 y<0O^P'O'l 

2MAbX J06400609-4327591 


ATA '\7"\7' 

GALAXY 


0.0570 


y 


J0652. 1+7425 


103.044 


74.425 


3.29+0.19 


17.1 


Mrk 6 


Syl. 5 


0.0188 


y 


J0656. 1+3959 


104.027 


39.986 


o oo 1 f\ o/^ 

2.29+0.26 


8.7 


TT/~1/^ 0/^01 

UGO 3601 


Syl 


0.0172 


y 


J0718. 0+4405 


109.517 


A A r\o A 

44.084 


1.67+0.24 


7.1 


OT\ T A O'VT' T0^1 OOO/^O 1 A ><OP'0^1 

2MAbX J 07180060+4405271 


Syl 


0.0610 




J0742. 5+4947 


115.644 


49.793 


O / \ O 1 n O "1 

2.98+0.21 


14.4 


Mrk 79 


Syl. 2 


O rAO oo 

0.0222 


y 


J0800. 1+2322 


120.032 


23.370 


1.62+0.24 


6.6 


SDSS J0759. 87+232448.3 


AT A ■\7" A r 

GALAXY 


0.0290 


y 


J0800. 3+2638 


1 OA nnn 
120.099 


26.648 


1.79+0.24 


/ .5 


IG 486 


byl 


0.0272 


y 


J0804.2+0506 


121.050 


5.101 


3.18+0.25 


12.6 


UGC 4203 


Sy2 


0.0135 


y 


J0811.1+7602 


122.798 


76.049 


1.15+0.19 


6.2 


PG 0804+761 


Syl 


0.1000 




J0814.4+0423 


123.600 


4.400 


1.26+0.24 


5.2 


CGCG 031-072 


Syl 


0.0331 




J0823. 2-0456 


125.800 


-4.947 


1.35+0.23 


6.0 


SWIFT J0823.4-0457 


Sy2 


0.0218 


y 


J0832.8+3706 


128.200 


37.100 


1.03+0.20 


5.2 


RB 0707 


Syl.2 


0.0919 




J0839.8-1214 


129.950 


-12.248 


1.28+0.21 


6.1 


30 206 


Syl 


0.1978 


y 


J0904.9+5537 


136.250 


55.632 


1.04+0.17 


6.0 


SWIFT J0904.3+5538 


Syl 


0.0370 


y 


J0911.5+4528 


137.898 


45.471 


1.23+0.18 


7.0 


SWIFT J0911. 2+4533 


GALAXY 


0.0268 


y 


J0918.4+1618 


139.615 


16.316 


1.65+0.21 


8.0 


Mrk 104 


Syl.5 


0.0292 


y 


J0921. 0-0803 


140.257 


-8.067 


2.59+0.20 


12.9 


SWIFT J0920.8-0805 


Sy2 


0.0198 


y 


J0923.8+2256 


140.962 


22.936 


2.05+0.20 


10.5 


MCG +04-22-042 


Syl.2 


0.0327 


y 


J0925.2+5217 


141.316 


52.285 


3.03+0.17 


18.1 


Mrk 110 


Syl 


0.0353 


y 


J0945.8-1419 


146.468 


-14.332 


1.28+0.21 


6.1 


NGC 2992 


Sy2 


0.0077 


y 


J0947.7-3056 


146.939 


-30.948 


11.50+0.22 


51.5 


ESO 434-40 


Sy2 


0.0085 


y 


J0947.9+0727 


147.000 


7.451 


1.12+0.21 


5.4 


30 227 


RG 


0.0860 


y 
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SWIFT NAME 


R.A. 

(J2000) 


Decl. 
(J2000) 


Flux 
flO-ii ess) 


S/N 


ID'' 


TvDe'' 


Redshift 


In BAT 22 months'^ 


J0959. 6-2250 


149.916 


-22.834 


4.44+0.22 


19.8 


NGC 3081 


Sy2 


0.0080 


y 


J1001.8+rx542 


150.453 


55.700 


1.43+0.16 


8.8 


NGC 3079 


Sy2 


0.0037 


y 


J1006. 0-2306 


151.500 


-23.100 


1.25+0.23 


5.5 


ESO 499-G 041 


Syl 


0.0127 




J1021. 7-0327 


155.450 


-3.450 


1.35+0.22 


6.3 


MCG+00-27-002 


Syl 


0.0409 




J1023.5+1951 


155.888 


19.864 


7.35+0.20 


36.8 


NGC 3227 


Syl.5 


0.0039 


y 


J1031.8-3451 


157.975 


-34.860 


4.71+0.26 


18.4 


NGC 3281 


Sy2 


0.0107 


y 


J1031.9-1417 


157.996 


-14.300 


2.13+0.23 


9.3 


H 1029-140 


Syl 


0.0860 




J1044.0-I-7023 


161.003 


70.400 


0.99+0.17 


5.9 


MCG+12-10-067 


Sy2 


0.0333 




J1046.5-I-2556 


161.649 


25.950 


1.12+0.19 


5.9 


UGC 05881 


GALAXY 


0.0200 




J1048. 5-2512 


162.149 


-25.200 


1.42+0.27 


5.3 


NGC 3393 


Sy2 


0.0125 




J1049.3-I-2256 


162.350 


22.950 


1.57+0.20 


8.0 


SWIFT J1049.4+2258 


Sy2 


0.0328 


y 


J1106.6-t-7234 


166.654 


72.571 


6.45+0.17 


38.0 


NGC 3516 


Syl.5 


0.0088 


y 


J1115.9+5426 


168.999 


54.450 


0.88+0.15 


5.7 


SDSS J111519.98+542316.6 


Sy2 


0.0703 




J1125.4-I-5421 


171.352 


54.351 


0.97+0.15 


6.3 


ARP 151 


Syl 


0.0210 




J1127.5+1908 


171.900 


19.148 


1.12+0.20 


5.6 


IRXS J112716.6+190914 


Syl 


0.1050 


y 


J1132. 7+5259 


173.188 


52.988 


1.01+0.15 


6.6 


UGC 6527 


Syl 


0.0277 


y 


J1136. 5+2132 


174.150 


21.548 


1.13+0.19 


5.9 


Mrk 739 


Syl 


0.0299 




J1139.0-3744 


174.764 


-37.741 


10.07+0.27 


37.9 


NGC 3783 


Syl 


0.0097 


y 


J1139.1+5912 


174.783 


59.212 


1.25+0.15 


8.1 


SBS 1136+594 


Syl.5 


0.0601 




J1139.4+3156 


174.869 


31.935 


1.00+0.17 


5.8 


NGC 3786 


Syl.8 


0.0089 




J1144.7+7939 


176.190 


79.662 


2.13+0.18 


11.9 


SWIFT J1143.7+7942 


Syl. 2 


0.0153 


y 


J1145.3+5859 


176.349 


59.000 


0.81+0.15 


5.3 


Ark 320 


GALAXY 


0.0099 




J1145.5-1825 


176.393 


-18.428 


2.84+0.27 


10.5 


2MASX J11454045-1827149 


Syl 


0.0329 


y 


J1148.9+2938 


177.230 


29.634 


1.06+0.18 


6.1 


MCG+05-28-032 


LINER 


0.0230 




J1158.0+5526 


179.502 


55.449 


1.04+0.15 


6.9 


NGC 3998 


Seyfert 


0.0036 


y 


J1201. 0+0647 


180.250 


6.800 


1.18+0.21 


5.6 


SWIFT J1200.8+0650 


GALAXY 


0.0360 


y 


J1203. 0+4432 


180.773 


44.534 


2.33+0.15 


15.1 


NGC 4051 


Syl.5 


0.0023 


y 


J1204.5+2018 


181.149 


20.301 


1.32+0.19 


7.1 


ARK 347 


Sy2 


0.0225 


y 


J1206.2+5242 


181.565 


52.710 


1.24+0.15 


8.3 


NGC 4102 


GALAXY 


0.0028 


y 


J1209.1+4700 


182.300 


47.000 


0.76+0.15 


5.0 


Mrk 198 


Sy2 


0.0246 


y 


J1209.4+4341 


182.370 


43.686 


1.56+0.15 


10.1 


NGC 4138 


Syl.9 


0.0030 


y 


J1210.5+3924 


182.633 


39.406 


24.60+0.16 


153.4 


NGC 4151 


Syl.5 


0.0033 


y 


J1210.6+3819 


182.667 


38.333 


0.96+0.16 


6.0 


LEDA 38759 


Syl 


0.0230 




J1217.2+0711 


184.300 


7.200 


1.21+0.20 


5.9 


NGC 4235 


Syl 


0.0080 


y 


J1218. 3+2950 


184.593 


29.839 


1.53+0.17 


9.1 


Mrk 766 


Syl.5 


0.0127 


y 


J1219. 0+4715 


184.750 


47.252 


0.96+0.15 


6.3 


NGC 4258 


LINER 


0.0015 


y 


J1222.0+7518 


185.503 


75.311 


1.27+0.17 


7.4 


Mrk 205 


Syl 


0.0700 


y 


J1225.7+1239 


186.447 


12.665 


12.58+0.19 


65.6 


NGC 4388 


Sy2 


0.0084 


y 


J1225.8+3330 


186.466 


33.513 


1.25+0.16 


7.7 


NGC 4395 


Syl 


0.0010 


y 


J1235.6-3955 


188.902 


-39.919 


10.21+0.26 


39.4 


NGC 4507 


Sy2 


0.0118 


y 


J1238.8-2718 


189.723 


-27.308 


4.39+0.28 


15.9 


ESO 506-027 


Sy2 


0.0240 


y 


J1239.0-1611 


189.769 


-16.196 


2.02+0.26 


7.7 


XSS J12389-1614 


Sy2 


0.0360 


y 


J1239.5-0520 


189.898 


-5.341 


4.52+0.23 


19.8 


NGC 4593 


Syl 


0.0090 


y 


J1246.6+5434 


191.661 


54.575 


1.34+0.15 


9.0 


NGC 4686 


GALAXY 


0.0168 


y 


J1302.8+1624 


195.700 


16.400 


0.90+0.17 


5.1 


Mrk 0783 


Syl.2 


0.0672 
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SWIFT NAME R.A. 


Decl. 


Flux 


S/N 


ID=- 


Type^ 


Redshift In BAT 22 months'^ 


(J2000) 


(J2000) 


(10-" cgs) 











J1306. 7-4024 


196.698 


-40.415 


2.37+0.27 


8.9 


ESO 32.3-077 


Syl.2 


0.0150 


y 


J1309. 1+1137 


197.279 


11.632 


2.19+0.18 


12.0 


SWIl^T ,11309.2+1139 


GALAXY 


0.0251 


y 


J1315. 4+4424 


198.852 


44.404 


1.28+0.15 


8.4 


IGR J13149+4422 


Seyfert 


0.0367 




J1322.3-1642 


200.591 


-16.716 


2.57+0.27 


9.5 


MCG -03-34-064 


Syl.8 


0.0165 


y 


J1325.4^4301 


201.366 


-43.017 


49.77+0.27 


187.6 


Cen A 


Sy2 


0.0018 


y 


J1334.8-2323 


203.700 


-23.400 


1.49+0.29 


5.1 


ESO 509-38 


Sy2 


0.0265 




J1335.7-3418 


203.944 


-34.302 


4.86+0.29 


16.6 


MCG -06-30-015 


Syl.2 


0.0077 


y 


J1338.1+0433 


204.547 


4.552 


3.63+0.20 


17.9 


NGC 5252 


Sy2 


0.0230 


y 


J1341.4+3022 


205.356 


30.369 


1.15+0.16 


7.2 


Mrk 268 


Sy2 


0.0404 




J1349. 5-3018 


207.390 


-30.304 


17.87+0.31 


58.2 


IC 4329A 


Syl 


0.0161 


y 


J1353.2+6919 


208.305 


69.327 


2.78+0.17 


16.7 


Mrk 279 


Syl.5 


0.0305 


y 


J1356.1+3835 


209.033 


38.583 


1.22+0.16 


7.7 


Mrk 464 


Syl 


0.0507 


y 


J1408.4^3024 


212.100 


-30.400 


1.65+0.32 


5.1 


PGC 050427 


Syl 


0.0235 




J1413.5-0312 


213.375 


-3.201 


14.38+0.24 


59.0 


NGC 5506 


Syl.9 


0.0062 


y 


J1418. 2+2507 


214.568 


25.133 


3.12+0.17 


18.2 


NGC 5548 


Syl.5 


0.0172 


y 


J1419. 5-2639 


214.893 


-26.663 


3.49+0.34 


10.4 


ESO 511-G030 


Syl 


0.0224 


y 


J1421. 6+4750 


215.420 


47.838 


1.03+0.16 


6.4 


QSO B1419+480 


Syl 


0.0720 


y 


J1424.3+2435 


216.100 


24.600 


0.89+0.17 


5.1 


NGC 5610 


GALAXY 


0.0169 




J1429.6+0117 


217.400 


1.300 


1.26+0.23 


5.4 


QSO B1426+015 


Syl 


0.0860 




J1436.5+5847 


219.149 


58.798 


1.37+0.16 


8.3 


QSO J 1436+5847 


Syl 


0.0312 


y 


J1441.2+5330 


220.300 


53.500 


0.85+0.16 


5.1 


Mrk 477 


Sy2 


0.0380 




J1442. 6-1713 


220.664 


-17.223 


4.83+0.34 


14.4 


NGC 5728 


Sy2 


0.0095 


y 


J1453.1+2556 


223.282 


25.936 


1.29+0.18 


7.0 


RX J1453.1+2554 


Syl 


0.0465 




J1504.2+1025 


226.073 


10.417 


1.51+0.22 


6.8 


Mrk 841 


Syl 


0.0364 


y 


J1515.4+4201 


228.868 


42.033 


1.05+0.18 


5.9 


NGC 5899 


Sy2 


0.0085 


y 


J1536.2+5753 


234.061 


57.890 


1.43+0.18 


7.9 


Mrk 290 


Syl 


0.0296 


y 


J1548. 4-1344 


237.106 


-13.749 


2.91+0.39 


7.4 


NGC 5995 


Sy2 


0.0251 


y 


J1554.8+3242 


238.700 


32.700 


1.04+0.20 


5.2 


2MASX J15541741+3238381 


Syl 


0.0483 




J1618.4+3224 


244.600 


32.400 


1.07+0.21 


5.1 


3C 332 


RG 


0.1500 




J1628.3+5147 


247.082 


51.793 


2.45+0.20 


12.4 


SWIFT J1628.1+5145 


Syl.9 


0.0547 


y 


J1653.2+0224 


253.319 


2.404 


4.20+0.35 


12.1 


NGC 6240 


Sy2 


0.0245 


y 


J1822.1+6421 


275.541 


64.361 


1.10+0.21 


5.3 


QSO B1821+643 


Syl 


0.2970 


y 


J1824.2-5620 


276.057 


-56.348 


1.98+0.29 


6.9 


IC 4709 


Sy2 


0.0169 


y 


J1835. 1+3240 


278.791 


32.683 


4.67+0.21 


22.0 


3C 382 


Syl 


0.0579 


y 


J1837.1-5922 


279.284 


-59.368 


1.79+0.28 


6.3 


FAIRALL 49 


Sy2 


0.0200 


y 


J1838. 6-6523 


279.658 


-65.394 


6.23+0.28 


22.4 


ESO 103-035 


Sy2 


0.0133 


y 


J1842.4+7946 


280.616 


79.771 


5.82+0.19 


29.9 


3C 390.3 


Syl 


0.0561 


y 


J1845.1-6223 


281.297 


-62.399 


2.52+0.28 


9.0 


ESO 140-43 


Syl 


0.0141 


y 


J1857.3-7827 


284.341 


-78.464 


1.86+0.26 


7.1 


LEDA 140831 


Syl 


0.0420 


y 


J1921.2-5840 


290.323 


-58.677 


3.54+0.28 


12.6 


ESO 141-55 


Syl 


0.0366 


y 


J1942. 7-1018 


295.680 


-10.316 


4.30+0.30 


14.3 


NGC 6814 


Syl 


0.0052 


y 


J2009. 1-6103 


302.289 


-61.064 


3.03+0.26 


11.5 


SWIFT J2009.0-6103 


Syl 


0.0149 


y 


J2018.3-5538 


304.598 


-55.649 


1.65+0.27 


6.1 


PKS 2014-55 


RG 


0.0600 




J2042.7+7508 


310.685 


75.136 


3.23+0.19 


16.9 


4C +74.26 


RG 


0.1040 


y 


J2044. 1-1043 


311.039 


-10.731 


5.62+0.29 


19.7 


Mrk 509 


Syl.2 


0.0344 


y 
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SWIFT NAME 


R.A. 

(J2000) 


Decl. 
{J2000) 


Flux 
(10~l^ cgs) 


S/N 




Type'' 


Redshift 


In BAT 22 months'^ ? 




Q1 Q 017 
OJ-O. Ul 1 


-o t .uuo 




J-O.O 






n ni 1 


y 


791 HQ 1 OQ'^Q 


O -1- ( . ouu 


-9.652 


1 f^fi-i-n 97 


5.9 


IN 91 07 nQ7 


T,TNP,R 






J2132. 1-3343 


323.028 


-33.727 


2.92±0.27 


10.7 


CTS 109 


Syl 


0.0297 


y 


J2136.0-6223 


324.006 


-62.400 


2.42±0.23 


10.6 


QSO J2136-6224 


Syl 


0.0589 


y 


J2138.8+3206 


324.713 


32.115 


1.29±0.20 


6.3 


LEDA 67084 


Syl 


0.0250 




J2200.7+1033 


330.199 


10.565 


1.76±0.21 


8.6 


SWIFT J2200.9+1032 


Syl.9 


0.0266 


y 


J2202. 1-3152 


330.526 


-31.878 


8.01±0.25 


31.8 


NGC 7172 


Sy2 


0.0087 


y 


J2204. 5-1-0335 


331.149 


3.600 


1.33±0.21 


6.3 


IRAS 22017+0319 


Sy2 


0.0610 




J2209. 5-4709 


332.387 


-47.166 


3.02±0.23 


13.4 


NGC 7213 


Syl.5 


0.0277 


y 


J2223.8-0207 


335.962 


-2.121 


1.93±0.22 


8.9 


3C 445 


Syl 


0.0564 


y 


J2235.8-2603 


338.966 


-26.054 


2.76±0.24 


11.6 


NGC 7314 


Syl.9 


0.0048 


y 


J2236.1-I-3357 


339.040 


33.952 


1.66±0.19 


8.8 


Arp 319 


Sy2 


0.0225 


y 


J2236. 8-1235 


339.223 


-12.599 


1.43±0.23 


6.2 


Mrk 915 


Syl 


0.0240 


y 


J2245. 7-1-3941 


341.449 


39.695 


1.71±0.18 


9.2 


3C 452 


Sy2 


0.0811 


y 


J2254. 1-1734 


343.535 


-17.578 


5.67±0.23 


24.7 


MR 2251-178 


Syl 


0.0640 


y 


J2258. 9-1-4053 


344.749 


40.899 


1.31±0.18 


7.2 


UGC 12282 


Syl 


0.0171 


y 


J2259. 54-2455 


344.899 


24.929 


1.51±0.19 


8.0 


LEDA 70195 


Syl 


0.0338 


y 


J2303. 2-1-0853 


345.809 


8.885 


3.87±0.20 


19.4 


NGC 7469 


Syl. 2 


0.0163 


y 


J2304. 7-0841 


346.194 


-8.686 


6.09±0.22 


27.9 


Mrk 926 


Syl.5 


0.0469 


y 


J2304.7-I-1217 


346.200 


12.300 


1.11±0.20 


5.6 


NGC 7479 


Sy2 


0.0079 




J2318.4-4221 


349.614 


-42.360 


4.09±0.20 


20.7 


NGC 7582 


Sy2 


0.0053 


y 


J2319.0-I-0014 


349.762 


0.241 


2.82±0.21 


13.5 


NGC 7603 


Syl 


0.0293 


y 


J2326.3-I-2154 


351.600 


21.900 


0.97±0.19 


5.1 


RBS 20005 


Syl 


0.1200 




J2342.0-I-3035 


355.500 


30.600 


1.08±0.19 


5.7 


UGC 12741 


GALAXY 


0.0174 


y 


J2358. 7-6052 


359.699 


-60.876 


1.10±0.19 


5.9 


PKS 2356-61 


RG 


0.0963 





''Sources with a SWIFT name were identified in the works of lTueller et al.l ||2008| . I2OO9I) . 

^R G are radio galaxies, while sources identified as Galaxies are candidate radio-quiet AGNs which are of the XBONG type (see lTueller et al] 
l2009h or for which an optical spectrum is not yet available. 

■^Is the source detected in the 22 months BAT survey of lTueller et al.l ||2009|') ? 



